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ABSTRACT 


Class-E amplifier operations are based on zero voltage, zero current 
switch turn-on. Forced turn-off of the switch is required, but the build-up 
of voltage across the switch is delayed. At “Optimum Operation” and 
under damped “Suboptimum Operation”, its efficiency is ideally 100%. 
Analysis of basic class-E amplifier have been done by (i)Using state equa- 
tions, (ii)Assuming constant input current and sinusoidal output current & 
(iii)Assuming constant input current & non-sinusoidal output current. The 
second analysis gives approximate design values of the amplifier and the last 
analysis gives more accurate design of basic class-E amplifier. The push- 
pull class-E amplifier have been analyzed by two methods, (i)Using state 
equations and (ii) Assuming constant input current and sinusoidal output 
current. The design equations obtained by the second method of analysis 
gives sufficiently accurate design values. 
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Chapter 1 


Introduction 


Class-E tuned amplifiers are high efficiency amplifiers. The high efficiency 
is achieved by reducing collector power dissipation of switching transistor 
used in the amplifier. The reduction in switching loss, reduces the heat sink 
volume and weight considerably or reduces the junction temperature rise, 
thereby decreases the transistor failure rate. 

At optimum operation of class-E amplifier, its efficiency is ideally 100%. 
The switching transition “off to on” is with the condition that collector-to- 
emmiter voltage and slope of the voltage waveform must be zero. “On to 
off” transition is with delayed voltage rise across the switch. Hence, for 
optimum operation, ideally there are no switching losses. If the switching 
• “off to on” transition time is significantly large, the power dissipation will 
be minimum because of zero slope of voltage across the switch. The output 
power, normalized to the product of switch peak voltage and peak current 
will be maximum at a switch duty-cycle 50%. In the thesis, the design 
of basic class-E amplifier and push-pull class-E amplifier have been done 
with 50% duty-cycle of operation. Basic class-E amplifier has very simple 
configuration, because of single switch. This amplifier output response has 
harmonics in which first two harmonics are significant. It requires harmonic 
filter for radio-frequency applications. The push-pull configuration of this 
amplifier have output response with no even harmonic components. 


Organization of Thesis 

Chapter-2 of the thesis is mainly dedicated to analysis and design of basic 
class-E amplifier. Section 2.1 describes the basic criteria imposed on switch 
voltage, circuit configuration of basic class-E amplifier and its operation, 
other possible circuit configurations of class-E amplifier and different appli- 
cation of class-E amplifier. Section 2.2 is describes exact analysis of basic 
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ciass-E amplifier using state equations. State equations of the amplifier for 
two modes are written and method of solution has been described. In sec- 
tion 2.3 based on solutions of exact analysis of section 2.2, simulation studies 
have been done. Curves describing for transient behaviour and steady-state 
behaviour of the amplifier response were plotted. Latter a study has been 
carried out on following: 


• Effect of change in power with supply voltage. 

• Effect of change in load resistance. 

• Effect of simultaneous change in load resistance and series tuned in- 
ductance. 

• Effect of change in power with frequency. 

Section 2.4 gives harmonic analysis of load current and voltage across 
the switch. Section 2.5 describes analysis of the class-E amplifier assuming 
constant input current and non-sinusoidal output current. First the analysis 
has been done by assuming output current to be fundamental. Latter, finite 
number of harmonics have been included in the analysis. A comparison of 
the two analysis with the exact analysis (described in section 2.2) is given. 
Section 2.6 gives design equations, which are derived from the analysis 
described in section 2.5. Approximate design equations were given based 
on the analysis, assuming sinusoidal output current. 

Modified procedure for more accurate design has been described, it is 
based on the analysis, assuming non-sinusoidal output current. The design 
method is for the following two options: 

1. Constant P 0 , Qi and f. 

2. Constant R, Qi and f. 

Chapter-3 deals with analysis and design of class-E push-pull amplifier. 
First section of the chapter introduces with the configuration and possible 
modes of operation of the class-E push-pull amplifier. Section 3.2 gives 
exact analysis of push-pull class-E amplifier with two modes of operation, 
using state equations. State equations for the two modes of the circuit 
operation were written and solved with a set of given circuit elements. 
Steady-state curves of the circuit parameters have been plotted. In section 
3.3, harmonic analysis of the output load current were carried out for steady- 
state and optimum operation. Section 3.4 gives the analysis of push-pull 
class-E amplifier, assuming constant input current and sinusoidal output 
current with two modes of operation. The results of the analysis of section 
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3.4 and section 3.2 were compared. In section 3.5 design equations for 
optimum operation of push-pull class-E amplifier were derived from the 
analysis given in section 3.4. Design example is given to design the circuit 
elements of the amplifier. 

Section 3.6 gives exact analysis of push-pull class-E amplifier with three 
modes of operation, using state equations. First order differential equations 
were written for the “three modes of operation” of the amplifier. These 
equations have been solved and steady-state circuit responses were plotted. 
Section 3.7 gives the harmonic analysis of output load current for three 
modes of operation. In section 3.8, analysis assuming constant input current 
and sinusoidal output current with three modes of 'operation is given. A 
comparison of the analysis of section 3.6 and the analysis of section 3.8 is 
also given. Section 3.9 deals with the design of class-E push pull amplifier 
for optimum operation at D=|. The design equations are derived from the 
optimum analysis of push-pull class-E amplifier for D=|. A design example 
is given to obtain circuit element values for optimum operation. 



Chapter 2 

Analysis & Design of Basic 
Class-E Tuned Amplifier 


2.1 Introduction 


2.1.1 Basic Criteria 

The class-E power amplifier operation is a circuit topology, which consists 
of a switch and a load network, with three basic criteria imposed on voltage 
across the switch. 

1. At the time of switch turn “on”, the voltage across the switch should 
be zero. 

2. The rate of change of voltage across the switch should also be zero, 
at the time switch turn “on” . 

3. At the time of switch turn “off’ , the voltage build-up across the switch 
should be delayed. 


2.1.2 Basic Class-E Amplifier and Its Operation 

The basic class-E amplifier circuit consists of a switch, a diode, a shunt 
capacitor, a series tuned circuit across a voltage source in series with a large 
inductor L x (RF Choke), as shown in Fig. 2.1. A switch driving signal has 
been assumed, which operates the switch, at desired operating frequency 
and duty-cycle. The capacitor C x , includes any capacitance inherent in the 
switch Si and diode D x , but its value is assumed independent of the voltage 
across it. Series inductor L 0 and capacitor C 0 may absorb any equivalent 
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inductance and capacitance of the load circuit respectively. The resistance 
R is sum of the, (1) input resistance of the load, (2) series a.c. resistances 
of L 0 and C Q , (3) equivalent lumped resistance of the a.c. losses in shunt 
capacitor C\, switch Si and diode D x . The large inductor L x in series with 
voltage source provides approximately a constant D.C. current. 



Figure 2.1: Basic class-E amplifier. 

The switch driving signal operates the switch S x periodically. There are 
two modes of the circuit, one during switch “on” , other during switch “off” 
time. At switch “on” mode; L 0 , C 0 and R forms a series resonant circuit 
and current source is shorted. During switch “off” mode; C x , L 0 , C 0 and R 
forms a series resonant circuit and current source forces current through the 
network. The diode D x does not allow a negative voltage across the shunt 
capacitor C\. A typical wave form for the two modes of the “optimum 
class-E operation” is shown in Fig. 2.2. 



Figure 2.2: A typical waveform of optimum class-E amplifier. 
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2.1.3 Different Circuit Configurations of Class-E Am- 
plifier 

1. Class-E tuned power amplifier with shunt capacitor and series tuned 
output circuit is shown in Fig. 2.1. This circuit will be analyzed in 
section 2.2. 

2. The push-pull configuration of the class-E amplifier is combination of 
two single-ended class-E amplifiers as shown in Fig. 2.3. This con- 
figuration provides symmetrical output waveforms and larger power 
output. The detail analysis of this configuration will be given in 
chapter-3. 



Figure 2.3: Push-pull class-E amplifier. 

3. Class-E tuned power amplifier with only one inductor L and one ca- 
pacitor C in load network is shown in Fig. 2.4. The resistor R is the 
load in which a.c. power is delivered. Principle of amplifeir operation 
and exact analysis is given in Ref. [4]. 

4. Class-E tuned amplifier configuration with finite D.C. feed induc- 
tance and no capacitance in parallel with switch is shown in Fig. 
2.5. This configuration is suitable for low frequency applications 
(10 KHz to 500 KHz). L\ is not an RF choke, but a finite inductance. 
The switch “off” operation is at zero magnitude and zero slope of 
switch current. But at the time of turn “on”, finite voltage exists 
across the switch. Thyristor can also be used, instead of switching 
transistor. The detail analysis of this configuration has been given in 
Ref. [7]. 

5. Dual configuration of basic class-E amplifier is shown in Fig. 2.6. 
The series inductor L includes the effect of transformer. All circuit 
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elements and responses are dual of the circuit shown in Fig. 2.1. The 
switch must now be operated at zero current. In Fig. 2.1, the net 
effect of series-tuned circuit at operating frequency is inductive, in 
Fig. 2.6. the net effect of parallel tuned circuit at operating frequency 
will be capacitive. This circuit configuration is discussed in Ref. [2]. 



Figure 2.4: Class-E circuit with one inductor and one capacitor in load 
network. 



Figure 2.5: Class-E amplifier with finite DC feed inductance and no capac- 
itance in parallel with switch. 



Figure 2.6: Dual of series tuned class-E amplifier with shunt capacitor. 

6. Parallel tuned class-E amplifier with series capacitor is shown in Fig. 2.7. 
The analysis of this circuit is given in Ref. [2]. Waveforms of this cir- 
cuit is given in Ref. [15]. 
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Figure 2.7: Parallel tuned class-E amplifier with series capacitor. 

2.1.4 Application of Class-E Amplifier 

Class-E circuit has been used in following applications. 

1. High efficiency solid-state VHF-UHF tuned power amplifiers Ref. [3,13]. 

2. High frequency electric process heating Ref. [19,20,21,22]. 

3. Electronic lamp-ballast Ref. [27]. 

4. Zero voltage switching or zero current switching class-E rectifier Ref. 
[14,15,16,17,18], 

5. Resonant regulated DC/DC power converters Ref. [23,24,25,26]. 

2.2 Analysis and Simulation of Basic Class- 
E Amplifier Using State Equations 

2.2.1 Analysis 

The basic class-E power amplifier to be analyzed is shown in Fig. 2.8. The 

analysis has been done with following assumptions, 

1. The active device Si acts as an ideal switch i.e., zero “on” resistance, 
infinite “off” resistance, and zero switching time. 

2. Diode has zero forward resistance, infinite reverse resistance and zero 
transition time from conduction state to blocking state or vice-versa. 

3. Inductor Li (RF choke) is large enough to provide constant input 
current at steady-state operation. 


4. All passive elements are ideal. 
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5. The shunt capacitance C\ is independent of switch voltage and in- 
cludes any capacitance inherent in switch and diode. 



Figure 2.8: The class-E amplifier with diode. 

The switch Si in Fig. 2.8 is operated periodically. The switch “on” 
interval is DT and “off” interval is (l-D)T. Where, D is switch duty-cycle 
and T is time-period of the switch operation. Hence, the circuit operates 
in following two modes. 


Mode-1 

Fig. 2.9 shows the equivalent circuit when switch S\ is “on” or diode D\ is 
“on” , since the circuit includes the diode across the switch to allow reverse 
•current. The current I s i, which equals (. I\ — 1 0 ), is carried by switch Si if it 
is positive and by D\ if it is negative. We have three first order simultaneous 
differential equations and voltage V cl is zero for this mode. 



Figure 2.9: Equivalent circuit of the class-E amplifier for mode-1. 


dIi = K 

dt L\ 

dip _ — Vco _ Rio 
dt Ido f . q 


( 2 . 1 ) 


( 2 . 2 ) 




10 


dVco Io 

dt C 0 


(2.3) 


Mode- 2 

Fig. 2.10 shows, the equivalent circuit for this mode when switch is “off” 
and diode D\ is reverse biased due to charge on C\. For this state, we have 
four first order differential equations. 



Isl| 


+ Vs 


I 11 

4 Io 

> 

Lo 


Icl 


* 


Vcl 





VCO — ! 



Figure 2.10: Equivalent circuit of the class-E amplifier for mode-2. 



dli _V t - Vd 

dt Li 

(2.4) 

dl 0 - 

__ V c i — V co RI 0 

(2.5) 

dt 

L 0 Lq 


C+. oN 

o 

II 

(2.6) 


dVd h - I 0 

(2.7) 


dt Ci 


Eq. (2.4), (2.5) and (2.6) are same as Eq. (2.1), (2.2) and (2.3) re- 
spectively with V c i = 0.0. Therefore, only one set of equations need be 
programmed for simulation. If during switch “off” interval, the voltage 
Vcl goes negative, the diode Di will conduct and circuit will revert to 
mode (1). 


2.2.2 Circuit Parameters 

Following circuit parameters have been chosen to simulate the circuit: 
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V s = 100V, Li = 5 mH, f = 50 KHz, D = 0.5, R = 3.089ft, C 1 = 0.203 fxF, 

C 0 = 0.13/zF, L 0 = 90.13 pH. 

2.2.3 Method of Solution 

Method of solution for obtaining steady-state characteristics by simulation 
is as following. 

1. It is assumed that, the circuit starts with switch ”on” at time t=0, 
the three first order simultaneous differential equations (2.1) — (2.3) 
are solved for interval between zero and DT with all initial conditions 
equal to zero. 

2. Solutions of I Q , V co and V c i at time DT are taken as initial con- 
ditions to start solution for switch “off” interval from DT to T. For 
this interval four simultaneous differential equations (2.4) — (2.7) are 
solved. 

3. Next cycle repeats with solutions of Ii, I 0 , V co and V& at the end of 
previous cycle, taking as initial conditions for the interval zero and 
DT. For interval DT and T step (2) is repeated. 

The above solution steps (2) and (3) repeat till steady-state is achieved i.e., 
•values of J l5 I 0 , V co and V c \ at time (n+l)DT and (n+l)T almost equals 
the respective values at nDT and nT. Where, n is the number of cycle of 
switch operations. 

Flowchart and computer program used to simulate the class-E amplifier 
is given in appendix (A). 


2.3 Simulation Study 

2.3.1 Transient Solutions 

Fig. 2.11 shows, how V cl and I 0 reaches from transient state to steady- 
state. With the circuit parameters given in section 2.2.2. It is shown that 
the steady-state is reached in around 400 cycles of operation. The plot 
shown is the curve for every tenth cycle. The circuit starts with switch-on 
for half time-period, then switch-off for next half time-period. At steady 
-state the peak value of V c i is 358.75V and peaks of I a are +34.92A and 
-3 1.38 A. For optimum circuit values, maximum stress on the switch is the 
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peak value of V ci . To start the circuit the supply voltage is applied at the 
same instant or after the switching signal starts. 

Fig. 2.12 shows V& vs I 0 curve at steady-state for one time-period. This 
figure is the outermost curve of Fig. 2.11. The switch turn-on occurs at t=0 
sec., at this instant V c \ and its slope is zero. The output current equals to 
source current value. The I 0 goes to its negative peak value and returns to 
a value close to source current in magnitude, but negative, at this instant 
switch is turned-off. The current I 0 and voltage V& starts increasing. The 
capacitor voltage V c i reaches to its peak value when current I 0 equals the 
source current value I\. Further, V cl starts decreasing and I 0 increases to 
a positive peak value. Now, both V& and I 0 decreases, till the switch is 
turned-on. At switch turn-on, current I 0 equals the source current I\ . 


2.3.2 Steady-state Solutions 

Taking the optimum values of the design components steady-state simula- 
tion results have been obtained. The curves of the steady-state solution 
shown in different plots have, the first half of the time-period as the switch 
“on” interval i.e., zero to DT and the second half of the time-period as 
switch “off” interval i.e., DT to (l-D)T. 

Fig. 2.13 shows curve of voltage across the switch (V cl ) for one time- 
period at steady-state. It has been observed that at the time of switch “off” 
the shunt capacitor voltage V& starts rising parabolically up to 358.75V, 
and at the time of switch “on” the voltage V c i decreases to zero, also the 
slope of V c i becomes zero. Here, the advantage of zero voltage, zero slope 
turn “on” is that even if switching-on time is a significant fraction of switch 
time-period, the switching losses will be minimum. 

Fig. 2.14 shows the curve of voltage across series capacitor C 0 for one 
time period at steady-state. The maximum voltage is 905V and minimum 
voltage is -723V. The two peak values are unequal because the circuit con- 
figuration for the two halfs of time-period are different. The average voltage 
value of this waveform is equal to supply voltage V s . 

Fig. 2.15 shows plots of source current Ii and the current which flows 
through either capacitor C\ or switch Si for one time period. For optimum 
operation at steady-state the current through diode is zero. The source cur- 
rent Ii is constant (17.35A approx) at steady-state, its magnitude depends 
on circuit parameters and magnitude of supply voltage source. At the time 
of switch turn “off” the switch current is slightly higher than twice of the 
source current (37.5 A). This requires forced turn “off” of the switch. But 
at turn “on” the switch current starts rising from zero value. The peak 
value of the switch current I s i pk is 48. 8A. 
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Fig. 2.16 shows output load current for one time- period. It is not pure 
sinusoidal. The two peaks of the current are not equal in magnitude. The 
two peak values are 34.92A and -31.38A. 
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Figure 2.15: Constant source current (II) and the current (Isl) which flows 
through the switch or capacitor. 



0 le-05 2e-05 


Time-period (sec) 

Figure 2.16: Output current for one cycle. 


2.3.3 Study the Effect of Change in Power with Volt- 
age Source 

Fig. 2.17 shows the variation of power with change in supply voltage mag- 
nitude. Power follows the square law with source voltage V s . Hence, the 
source voltage V s and output current I 0 has linear relation. 
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Figure 2.17: Supply voltage vs power output in load resistor. 

2.3.4 Effect of Change in Load Resistance 

Fig. 2.18 is plot of load resistance vs power output in the load. If the 
resistance R is decreased from its optimum value 3.089f2 to a value very near 
to zero, the power output in this will also decrease. If it is increased beyond 
the optimum value the power output decreases because shunt capacitor 
discharges at the time of turn “on” . 



Figure 2.18: Curve for change in output power with output resistance. 

Fig. 2.19 Shows shunt capacitor voltage Va for three resistance values. 
The plot is for one time-period at steady- state. If resistance R is decreased 
below the optimum value the circuit becomes underdamped, i.e., V cl goes 
negative at the time of switch turn “on” , but the diode D x starts conduction 
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and holds V& at zero (plot (a) with R=0.l£2). The increase in resistance 
R beyond optimum value the circuit becomes overdamped, i.e., at the time 
of switch turn “on” V cl is at some positive voltage. At the time of switch 
rum “on” V cl discharges through switch, resulting in power loss. 
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Figure 2.19: Variation of voltage across shunt capacitor with change in load 
resistance from its optimum value. 

During optimum operation of class-E amplifier, the diode Di does not 
conduct at steady-state. Because at time of switch turn “on” the V c \ and 
its slope both are zero i.e., V& never goes to negative value. The efficiency 
of the circuit at optimum operation is nearly 100 percent. 

During the suboptimum operation of the class-E amplifier, if V c \ is at 
some positive value at the time of switch turn “on”, the diode never con- 
ducts during steady-state. But this forces the discharge of capacitor C\ 
through switch S\ at turn “on”. This results in loss in energy and it may 
cause very large current through the switch, which can damage it. Since 
the circuit has ideal parameters, the loss which affects the efficiency is only 
the discharge of capacitor at the time of turn “on” . 

Suboptimum operation of class-E amplifier, in which voltage V& tends to 
go to negative, forces diode D\ to conduct during steady-state. The switch 
is turned “on” at zero voltage, and switch current will start building up 
after the conduction of diode is over. The frequency, duty-cycle and load 
resistance variation is possible independently up to some extent without 
affecting the efficiency. The efficiency of the circuit is nearly 100 percent. 
Because the switch-on takes place at zero voltage and non-zero slope (finite 
turn-on current). 

Fig. 2.20 is the plot of resistance vs peak voltage across shunt capacitor 
C\. This plot helps in selecting the device voltage rating, where the load 
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resistance value changes (this is usual in induction heating loads). Here, the 
Maximum voltage value is peak of any voltage that occurcd from transient 
stale to steady-state. The steady-state peak voltage value may be less than 
the peak voltage (lining transition from initial state to steady-state. For 
example. Fig. 2.1‘) shows peak I)., - M2 V for R. = 0.1 12, and Fig. 2.20 
and 2.21 shows maximum L ;I 700 1’. 



It (Olmi) 


Figm e 2.20: IVnk voltage across switch during transients with different 
values of load lesistauce. 


Fig. *’.‘’1 shows. the* plot of peak shunt capacitor voltage Ki,,,,,, vs respec- 
tive switching cycle. The plot is for lirst 1000 cycles. The circuit operation 
reaches steady-state in around .‘100 to700 cycles. It is observed that if the 
load resistance is decreased below a certain value (1.0 12 for this circuit), the 
shunt capacitor ‘voltage become underdamped. 'Phis puts very largo voltage 
stress across the switch. 'Phis is taken care in selecting switch ratings. 


2.3.5 Change in Load Resistance and Series Tuned 
Inductance 

When the chuss-E amplifier is applied for power supply of high frequency 
induct ion heating, the change in load resistance R and change in load induc- 
tor l.„ is very common. Here, a study has been done for this case. The load 
resistance, is the resistance of induction heating “job”. The inductance L 0 
includes I lie inductance of “job work”. The resistance R generally increases 
with temperature', and olfertive inductance generally decreases. Nature 
of change. of 11 and L„ depends on the material and its way of processing' 
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Figure 2.21: Peak voltage across switch during transients vs switch cycle 
for different load resistances. 


with temperature. In the present study it is assumed that increase in R 
and decrease in L„ is linear, as shown in Fig. 2.22. 

Fig. 2.22 shows the increase in resistance from 1S2 to 3 il and decrease 
in L„ from 102 fiH to 90 //if. The range of change in R and L n is chosen 
such that the switch turn “on” and turn “off” should be at zero voltage. 



Figure 2.22: Assumed characteristics of Lo vs R with temperature during 
induction heating. 

Change in power output is shown in Fig. 2.23 from operating point with 
R = 1 12, L„ = 102 fi.H to 7? = 3 S 2, L„ — 90 /iH. It is observed that the 
power output, to the load will increase as resistance increases. The desired 
power output can be adjusted by changing supply voltage. 

Fig. 2.24 shows evidenco(e.g.,for three cases) that voltage across switch 
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will be zero both at turn “on” and at turn “off” . This plot is for steady-state 
and one time-period. 

t'i / r \ i 


Power 

(Watts) 



9e-05 9.24e-05 9.48e-05 9.72e-05 

L 0 (Henry) 


9.96e-05 0.000102 


Figure 2.23: Change in output power due to change in Lo and R according 
to characteristic shown in Fig. 2.22. 


V c i 

(Volts) 



Figure 2.24: Voltage across the switch at steady-state with variation in R 
and Lo. 

Fig. 2.25 provides help in selecting peak voltage rating of the device. 
This plot is for any peak voltage occured during transient state to steady- 
state transition. Although the steady-state peak voltage will be lower than 
the peak voltage during transient, the device should be capable of with- 
standing any peak voltage occuring across it at any time. Here, the vari- 
ation of maximum voltage across the switch is plotted for transition from 
operating point with R = 1 Q, L 0 = 102 jiH to R = 3 Q, , L 0 = 90/jlH. 
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(Volts) 
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9e-05 9.24e-05 9.48e-05 9.72e-05 9.96e-05 0.000102 

L 0 (Henry) 


Figure 2.25: Peak transient voltage across the switch. 


2.3.6 Change in Power with Change in Frequency 


If all circuit element values are held constant, and resistance R is chosen 
1.0 f2, there is a range of frequency (48.5 KHz to 53.7 KHz) in which it 
can be varied without affecting efficiency of the tuned power amplifier. The 
small change in frequency results large change in power output as shown 
in Fig. 2.26. This can be applied in some induction heating loads, where 
change in resistance and change in inductance of work coil is negligible. 
This can also be used to regulate power output, if input supply voltage 
fluctuates. 


Power (Watts) 



48500 49000 49500 50000 50500 51000 51500 52000 52500 53000 53500 54000 


Frequency (Hz) 

Figure 2.26: Change in power with frequency. 
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2.4 Harmonic Analysis 


The class-E circuit is single ended and has two modes of operation. It 
may contain harmonics in its output waveform. To obtain harmonic com- 
ponents in output current I 0 , and shunt capacitor voltage V c i, we use dis- 
create fourier transform method. To obtain harmonic components, we use 
NAG subroutine C06EAF. This subroutine is based on Discreate Fourier 
Transform. 


2.4.1 Harmonic Analysis of Output Current 

The harmonic components of output current I 0 shown in Fig.2.16 for op- 
timum case have been obtained. Peak magnitudes and phase angles of 
different harmonics are shown in table-2.1. It is observed that first three 
harmonics are significant, and rest of the harmonics are less than one per- 
cent of the fundamental. Phase angle of fundamental is 145.55 degree and 
peak magnitude is 33.15A. The phase angles of even harmonics are slowly 
decreasing from 221.19 degree, as the order of harmonics increases. 


Table-2.1 

Harmonics and phases of Output Current 

Harmonics 

Peak-magnitude (Amp.) 

Phase angle (Degree) 

1 

33.146 

145.55 

2 

1.925 

221.19 

3 

0.294 

-1.82 

4 

0.126 

201.05 

5 

0.0562 

-0.916 

6 

0.0335 

194.33 

7 

0.0198 

-0.443 

8 

0.0136 

190.98 

9 

0.0092 

-0.124 

10 

0.0068 

189.01 


2.4.2 Harmonic Analysis of Shunt Capacitor Voltage 

The harmonic components of shunt capacitor voltage V c \ shown in Fig.2.13 
for optimum case have been obtained. Magnitudes and phase angles of dif- 
ferent harmonics are listed in table-2.2. The average of the shunt capacitor 
voltage is equal to the supply voltage. First three harmonics are significant, 
and rest are very small. 




24 


Table-2.2 

Harmonics and phases of Shunt Capacitor Voltage 

Harmonics 

Peak-magnitude (Volt.) 

Phase angle (Degree) 

0 

99.999 

0.000 

1 

82.119 

196.161 

2 

42.838 

-52.980 

3 

11.301 

85.659 

4 

6.748 

-70.876 

5 

3.843 

87.469 

6 

2.775 

-77.125 

7 

1.934 

88.215 

8 

1.521 

-80.318 

9 

1.164 

88.624 

10 

0.962 

-82.252 


The shunt capacitor voltage will always be equal to the voltage across 
the series tuned impedance Zo. Hence, the average D.C. voltage blocked 
by series capacitor Co will be equal to the supply voltage Vs. The different 
harmonic components of output current can also be obtained by divid- 
ing the different harmonic components of shunt capacitor voltage by series 
impedance Zoh. Where Zoh is for respective harmonic frequency. 


2.5 Analysis Assuming Constant Input Cur- 
rent and Non- Sinusoidal Output Cur- 
rent 


A simple equivalent circuit of class-E amplifier is shown in Fig. 2.27. Anal- 
ysis of the circuit has been done with following assumptions. 

1. The RF choke Zq in series with supply voltage source provides ap- 
proximately constant D.C. input current Ij. 

2. The output current I 0 {9 ) is non-sinusoidal at the operating frequency 
of the amplifier. But average of the output current I 0 {9 ) is zero, at 
steady-state, because of output capacitor C 0 . The quality-factor of 
the series tuned output circuit is finite. 

3. The switch and diode have ideal characteristics. 

4. All passive elements are ideal. 
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In the equivalent circuit shown in Fig. 2.27 the series reactance jX is 
produced by series connection of L 0 and C a (Fig. 2.1) at operating frequency 
of the tuned amplifier. If the operating frequency is equal to the resonant 
frequency of L 0 and C 0 the reactance jX equals zero. 



Figure 2.27: Equivalent circuit of class-E amplifier. 

Let the voltage across resistance R be V 0 (9), voltage across jX is V x [9) 
and voltage across the shunt capacitor is V c \ (9). The current through series 
tuned circuit is I 0 (6). Hence, current I $ i (9) (shown in Fig. 1.27) is equal 
to the difference between the input current I\ and output current I 0 (9). 

We start the analysis, first by assuming the output current I a pure 
sinusoidal, then we impose the conditions of optimum operating conditions 
to obtain the circuit parameters. Later a general analysis have been done 
by assuming harmonics in output current. 


2.5.1 Analysis Assuming Output Current to be Fun- 
damental Sine- Wave 

If the quality factor of the series tuned circuit is very high the second and 
higher order harmonics of output current will be very small w.r.t. to fun- 
damental component. The analysis assuming only fundamental component 
in output current is simple but approximate. 

Suppose, the sinusoidal current I 0 {9) is given by, 


Io{8) = I om sin [9 + (j>) (2.8) 

Where, I om is the peak-current value, <fi is the phase angle w.r.t. to some 
arbitrary reference, and 9{= cot) is angle in radian corresponding to oper- 
ating frequency of tuned amplifier. The instant when switch is turned “on” 
is taken as the reference for phase of current I 0 {9) (Fig. 2.2). 
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When switch Si is “off 5 , capacitor C\, equivalent reactance jX and 
resistance R are in series. V c i (9) is given by, 


v c i(d) = v x (e) + v 0 (e) 

Hence, from Eq. (2.8) and (2.9), we have 


(2,9) 


V a (0) = I m vOf 2 + R 2 sin (« + <)> + if,) (2.10) 

Where, ip = arctan (X/R). 

The voltage V c i{9) exists only during switch is “off” , and diode D\ is 
not conducting. The average voltage across the source inductor L x is zero, 
because it carries a constant D.C. current. Hence, the average voltage across 
the shunt capacitor Ci, will be equal to supply voltage. 

The circuit model of the class-E tuned amplifier, chosen for approxi- 
mate analysis contains all lossless circuit elements. The loss occurs only 
in condition if the switch turn “on” takes place at some finite capacitor 
voltage value. The class-E circuit will operate at 100 percent efficiency, if 
the switch voltage drops to zero at the instant of the switch turn “on” . 

Power output increases as duty-cycle is increased, but the “peak-power 
output capability” for a given switching device can be achieved by operating 
the class-E amplifier at 50-percent duty-cycle and at optimum operating 
conditions, Ref.[2], At optimum operating point voltage across the switch 
does not become negative, hence, the circuit will work without diode D\. 
For very high frequency operation, if speed of diode is slower than the the 
switch speed, the optimum operating point is chosen. 

For the analysis of class-E tuned amplifier at optimum operation, we 
will choose duty-cycle D=0.5 and zero magnitude and zero slope of switch 
voltage at switch turn “on”. 

The switch Si is “on” for time duration zero to DT (i.e., 0 to it) and 
it is “off” from DT to T (i.e., 7r to 2ir). V c i (9) is non zero during switch 
“off” time (Fig. 2.2), and current I s i will flow only through the capacitor, 
because diode Di will not conduct at optimum operation. Equivalent circuit 
of Fig. 2.10 applies to the switch “off” interval. 

Hence, I s i(9) is given by, 

Isi{9) = h-IomSm {9 + (p) (2.11) 


Also, voltage V c i (9) is given by, 
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Va (0) = J [h ~ lom sin (a + 4>))da (2.12) 

Where angle 9 can be any value between 7r to 27 r. The Eq. (2.12) can 
be written as, 


Wl (9) — - 7r) + lorn cos (6 + <f>) + I om COS(f)] (2.13) 

The voltage across inductor equal V s , when switch Si is closed and 
(V $ — V c i) when the switch is open. In order that the average inductor 
voltage be zero in steady-state, the average value of V& must equal V s . 
Therefore 


V s = ^~ rv c i(9)de (2.14) 

Z7T J 7r 

Using Eq. (2.13) and (2.14), we get the following result. 

Ehau — Tr[di7r “b sin (p -f- COS (p) (2.15) 

AnwUi 

Average value of V cl (6) is equal to V s , 

V s = - 1 - - [Ji7T 2 + AI om sin (j) -I- 2 ITlorn cos <f>] (2.16) 

AnwCi 

Using Eq. (2.13), we obtain fundamental frequency component of V c i(9), 
as below, 


V& fundamental (9) 


[{— Ii7r — 2 Iom cos 4> — 7r ^° m sin <p} sin 6 


TTwCi 


+{2Ii + — cos 4>} cos 9] 


(2.17) 


From Eq. (2.10), we have, 

V* ( 9 ) = IomZ 0 cos (cf) + ip) sin 9 + I om Z 0 sin (<p + ip) cos 6 (2.18) 

The fundamental frequency co mponent of V& (9) of Eq. (2.17) is equal 
to the voltage drop across Z 0 (= \fXP- + W). Hence, from Eq. (2.17) and 
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(2.18), we equate coefficients of sin# and cos#, respectively and obtain 
following two equations; 


(wC\X - -)I 0 m sin </>-(- + wCiR)I om cos 6 = h 
Z 7T 


(2.19) 


1 2 1 

wCiRIom sin<^> + (-- + wCiX)I om cos (f>= — 

2 71 


From Eq. (2.19) and (2.20), we have 


( 2 . 20 ) 


. Il [2a>£iR +w c 1 X + ±- \] 

I om sin 4, - (toCiX _ l )2 + wCl R{ 1 + wCl R) 

omCOScp ( wCiX _ly + w C lR (%+wC 1 R) 


( 2 . 21 ) 


( 2 . 22 ) 


From Eq. (2.21) and (2. 22), we have, 


4> = arctan 


[B+^ + 


TTwCl 4wC 1 j 


[V-f 


2 wCi - 


(2.23) 


Generally, value of (f) obtained by Eq. (2.23) will be in either in first 
quadrant or in fourth quadrant. If I om cos o is negative, tt radian is added 
to (j) to make it in second or third quadrant. 

If current source I x is realized by a very large inductor in series with 
voltage source Vs, the value of I x can be obtained from Eq. (2.16), because 
I om sin 4> and I om cos 4> are known. 

Once, I x , Iom, and (j) are known, the output current I 0 (#) and voltage 
V cl (#) can be found at any instant of time. 

■ From Eq. (2.13) at # = 2n, we get, the value of V c i(#) as, 


v cl / e=2n = + 2 I om cos <S>) (2.24) 

Consider, peak value of V c i (#) occurs at # = 9 vmax , from Eq.(2.13), we 
get, 


dVjfi) 

dd 


/e 


=e v 


— 0 — I\ Iom sm(9 vma x ffi b) 


(2.25) 
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Since, values of A, I om and © are known, we can solve the value of 6 vmax 
from Eq. (2.25). Hence, V c i max can be obtained from Eq. (2.13), using 
6 = d VT nax- 


Solution for Optimum Case 

At the time of switch turn £: on”, the V c i(9) should be zero. Hence, from 
Eq. (2.13) at turn “on”, 


Vci( 6 ) / 6 = 2 it — 0 — [/i(2?r — 7r) + I om cos (27 r + <fr) + I om cos <f>] (2.26) 


which gives, 


cos 4> - 


ttIi 


2/„ 


(2.27) 


The slope of V c i ( 9 ) at 9 
(2.13), 


: 27 r must also be zero. Hence, again from Eq. 


dV c i(0) 

d9 


/ 9=2-k 


= 0 = [Ii - I om Sin (27 r + <j>)] 


Eq. (2.28) gives, 


(2.28) 


sin (p = (2.29) 

J-om 

Eq. (2.27) and (2.29) gives, 

4> = 2.57468 c (2.30) 

From equation (2.16), (2.27) and (2.29), we obtain, 


I\ — irwCiVs 


The output power in the load resistor is given by, 

■V iLR 


(2.31) 


2 


(2.32) 
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The input power is simply, 

Pi = V S I X (2.33) 

Hence, efficiency is given by, 

(2M) 

The maximum possible switch current is given by, 

* 

I sl(pk ) — l"l t Iom (2.35) 

From equations (2.29), (2.30) and (2.35), we get, 

Isi( P k) = 2.862 h (2.36) 

From equation (2.13) , at peak value of V c i(9 ), 

djj V' e=e vma x — 0 = J"i — hm sin ( 9 vmax + <^>) (2.37) 

From equations (2.29), (2.30) and (2.37), we obtain, 

Ovmax = 4.27539 c (2.38) 

Substituting value of 9 vmax for 9 in Eq. (2.13) and using Eq. (2.29), 
(2.30) and (2.31), we have V c i(9) at 9 max as, 

Kimax = 3.562V, (2.39) 

2.5.2 Analysis Assuming Output Current to be Non- 
Sinusoidal 

Suppose, the output current I 0 (9) is non-sinusoidal and is represented by 
sum of its k-harmonic components (where, k is an integer and k — > co).as 

k 

m = I° m n sin (n9 + <j > n ) 


(2.40) 
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Where. Iom n is the peak current value, and <f> n is the phase angle of n th 
harmonic current, w.r.t. to some arbitrary reference, and 9(— wt) is angle 
in radian. The instant when switch is turned-on is taken as the reference 
for phase of current I a (9) (ref. Fig. 2.2). From Fig. 2.27, V cl (0) is given by, 

V cl (9) = V x (9) + V 0 (9) (2.41) 

The equivalent reactance for n th harmonic will be jX n . Hence, from Eq. 
(2.40) and (2.41), we have, 


V cl (9) = V cl „ + Iom n ^xf+B}sm(n9 + + $ n ) (2.42) 

n= 1 

where, t/) n = arctan (^ ) . 

The current I s \ (8) (ref. Fig. 2.27) is given by, 

k 

I Si (#) = h ~ X ) Iomn sin ( n61 + <f>n) (2.43) 

71=1 

Also, voltage V c \{9) is given by, 

V e i (9) = — Iom n sin{na + <j) n )]da (2.44) 

WL/ 1 71=1 

This voltage V c i (9) exists only when switch is off. From Eq. (2.44), we 
have, 


V cl {9) = -^[h(9 - tt) + Y, —{cos (n9 + cf> n ) - (-l) n cos (2.45) 
wUi n=1 n 

The average inductor voltage is zero at steady-state, hence, the average 
value of V c i will be equal to Vs. Therefore, 

1 7-277 . 

V s = ~ V cl (9)d9 (2.46) 

Z7T J'K 

Using Eq. (2.45) and (2.46), we have, 
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x r 1 r 71 " 2 r i V'' r f 1 - (-1) 71 • , 

Vc = — I— h + Iom n{ — sm (p n 


2 ttiuC: 1 2 


n= 1 


n 2 


(-l) n7r 

n 


cos <J>n}] (2.47) 


If the impedance of series tuned output circuit is Zo m , for m — harmonics, 
we can obtain V cl (9), using Eq. (2.40), 


k 

v cl (6) V co av ™f" ^ ' I OTTiyyiZi Ofji sin(m0 H” (pm 'Pm') (2.48) 

m—l 

For m th harmonic component, ip m = arctan (2V); and average voltage 
across series capacitor Co is V COav (where, Vco av = V s ). From Eq. (2.45), 
V c i (&) can be represented by sum of its Fourier Components, as 

k 

Vci (0) = Vcl av + [a m cos (md) + b m sin {m6)\ (2.49) 

m— 1 

Where, V c i av is equal to supply voltage Vs, and 


1 r {l-(-l) m } r tt _ 

&m — /~t l i A d „ Iorn m cos (p T 

■kwCi m 2 2m 

k _ i} 

+ > 5 5 Iom n sm <p n \ 

m l — n 2 


n=l,n^m 


(2.50) 


1 r TT _ 7T . , X m{{- l) m+n - 1} 

ixwcP m 1 2m ommSm m + n= ^ m n(m 2 - n 2 ) 

+ jz WILzi ~ 1)m h om„ cos K] 


n = 1 


mn 


Iom n cos (j) n 
(2.51) 


Eq. (2.48) can also be written as, 


k 

Ki (9) = V ctav + [{Zo m Iomm sin (<f> m + ip m )} cos (m9) 

m= 1 

+{Zo m Iom m cos (m<j> + Vm)} sin(m0)] (2.52) 

Where, V cXav is equal to Vs. Since, Eq. (2.49) and (2.52) represent the 
same waveform, hence, 
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a m = Zo m Iom m sin(0 m 

' 1pm) 

(2.53) 

bm Zo m I om m COS (0m 

+ Ipm) 

(2.54) 


In Eq. (2.40), output current I 0 (6), consists of k — harmonic compo- 
nents. Hence, there will be 2k unknowns i.e., k — peak currents Iom n , and 
k — phase angles 4> n . If we choose, k — harmonic components of V cl (6) (i.e., 
m=l to k), from Eq. (2.50), (2.51), (2.53) and (2.54), we will have 2k equa- 
tions. These equations can be written in the form of linear equations, with 
assuming (Iom n sin <j> n ) and (7 om n cos <f> n ) as unknowns. 

Using Eq. (2.50) and (2.53), the m th equation can be written as, 


k / 2 > ) m + n — 1 

('KwCiR)Iom m sin <p m — ]T 1 — 15 — 7om n sin © n 


m * — n i 

7T 1 - (-l) m 

+{'KwC\X m — - — )Iom m cos 4> m = r R (2.55) 

TTr 


n=l, 

JR, 
2m J 


Using Eq. (2.51) and (2.54), the (m 4- k) th equation can be written as, 


tt (-l) m {1 - (- l) m } 

(- , KwC\X rn )Iom rr i sin 6 m 4- WwCiR ■= -}Iom m cos <b r 

2m m 2 

± + (- ir j-Eg Uom ,co S ^ - (2.56) 




mn 


m 


Eq. (2.55) and (2.56) each forms (m=l to k) k-equations. Hence, there 
are 2k equations and 2k unknowns (i.e., Iom n sin <fi n and Iom n cos <p n for 
n=l to k). Hence, all peak current values Iom n and <fi n can be solved for 
given value of operating frequency, input current source value R and circuit 
elements. 

If current R is unknown, we will have total (2k+l) unknowns. Thus, 
we will have to include one more equation i.e., Eq. (2.47). Hence, from Eq. 
(2.55), (2.56) and (2.47); we can write following matrix equation: 


Amn -Bmn 


' R 


’ E m ' 

Ym C mn D mn 


I OTTljji COS (bm 

= 

Fm 

i 

£ 

i 


I OTTIjji sill (prn 


G 


(2.57) 
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Where, 



" Xi ' 


’ Vi ~ 

x m = 


II 



. x k . 


. y k . 






&1 


/ 1 

v n = 

Vi . 

. . v k 

; E m = 

. e k . 

- F = 

5 1 m — 

_ fk . 



a n 

&ln 

Qifc 

Amn — 

&ml • 1 

&77172 

• • • 


. Ski 

&kn 



' b u 

bln 

1 

i— 1 

r-O 

Bmn — 

bml 

• • b mn 

- • • b m k 


. &Jfcl 

£ 

bkk . 


Of\Yv ^1V\VV <X-r<_ V-X k. W\<CtVv'ca>f. 

Here, both m and n vary from 1 to k. where, k is number of harmonics 
considered. Different elemental values of matrix is given as below. 

For m=l to k, 



“ 2 > tr TrL m 

G = 2'kwCiVs, f m = 0.0 


For n=l to k, 


V, = 


i - (-1)" 

r?2 
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V 

v n 

&nvn — 


_ (-1 fz 


, /or n 7 ^ to. 


n 

{(— 1)^+” — ]L} 

m 2 — n 2 
7rwCiR, for n — m. 

0.0, /or n 7^ m. 

7 r 

7rwCiX m — - — for n = m. 
2m 

0.0, /or n^m. 

7 r 

irwCiXm , for n — m . 

2m 


d -v\ v* &mzi — 


m{(-l)”+"-l l-(-l) m , , 

L ( o o'! "f" ( •*-) J i J or n 

n{m 2 — n 2 ) mn 


m. 


1 _ 

= [TnoCi-R ~ (-l) m ], /or n = m. 

m 2 


For given values of supply voltage Vs, operating frequency f, circuit 
element values Ci, L 0 , C 0 and R, we can find out I\ and output current 
Io(8) (by solving the matrix Eq. (2.55) ). Once, output current is known 
shunt capacitor voltage V c i(0) can be obtained at anv time instant from Eq. 
(2.45). 

If the circuit is operating at optimum condition, the input power will 
be equal to the output power. For a given value of power output Po, input 
current Il(=Po/vs) can be assumed constant. 

The flow-chart and program of the analysis assuming non-sinusoidal 
output current, is given in appendix (B). 


2.5.3 Comparison of Analysis Results With Results 
of Simulation 

For a given value of circuit elements, we can get solutions for circuit re- 
sponses by the two methods of analysis, (1) analysis assuming output cur- 
rent sinusoidal, (2) analysis assuming non-sinusoidal output current. There 
are several other methods, using which the class-E circuit can be analysed 
exactly (e.g., time-domain analysis and s-domain analysis). The importance 
of harmonic analysis method is that it leads to an accurate design method 
for optimum operation of the class-E amplifier. We will take a set of cir- 
cuit element values, with which the circuit operates at optimum operating 
conditions. The responses of the circuit will be compared with results of 
simulation and the results obtained by the two analysis methods described 
in previous sections. The circuit values are given below: 



v s = 100V, Li = 5 mH, f = 50Affz, D = 0.5, A = 2.7484ft, 
Ci = 0.2292 ./xF, C 0 = 0.1317/xF, L 0 = 87.4843/Ur. 
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Fig. 2.28 shows the plot of three curves of V&. Curve 1 has been ob- 
tained by calculation of V cX (6) from 7r to 27 t, using Eq. (2.13). This equation 
has been derived assuming output current I 0 pure sinusoidal (i.e., assuming 
fundamental component only). This analysis gives approximate solution. 
The peak value of V ci is 343.76V, and at 6 = 2-k its magnitude is 6.3V 
with negative slope. Curve 2 has been plotted by calculating V cl (9) from 
7r to 27 r, using Eq. (2.45). For this case, the output current I a is assumed 
non-sinusoidal. While calculation, first ten harmonics were considered in 
the I 0 . Since, magnitudes of harmonics beyond third harmonics, are very 
small. Hence, assuming ten harmonics in the calculation gives almost same 
result, as it is obtained by exact analysis shown by curve 3. The peak value 
of V c i in curve 2 is 358.87V, and at 9 — 2tt its magnitude and slope is almost 
zero. Curve 3 and curve 2 are almost same. The peak value of V cX in curve 
3 is 359.45v, and at 9 = 27 f its magnitude and slope are almost zero. 
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Figure 2.28: Plots of V cX to compare with simulation result. 

Fig. 2.29 is plot of I\ and I 0 . The input current I x has been assumed 
20. 0A in the analysis, and the plotted I x has been obtained by simulation. 

■ Three curves of I 0 has been shown in the figure. Curve 1 has been obtained 
by calculating I a from Eq. (2.8) i.e., assuming pure sinusoidal output cur- 
rent I 0 . Its magnitude at 9 = 2-n is 21.83A, and peak values are -37.46A 
and +37.46A. Curve 2 is plot of output current which has been calculated 
from Eq. (2.40), and assuming first ten harmonics only. The analysis as- 
suming output current with first ten harmonics gives almost the same result 
as obtained by exact analysis. This analysis gives the value of I 0 19.46A at 
9 = 27r. Its peak values are -38.49A and +38.49A. Curve 3 shows the plot 
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of output current I 0 obtained by simulation result (i.e.. exact analysis of 
the circuit). At 6 — 2-rr the magnitude of I 0 is 19.86A. Its peak values are 
-36.33A and +40.02A. 

The values of I 0 at switch turn-off, given by curves 1, 2, and 3 are - 
22.31A, -19.98A and -22.63A. Current I sw \ (equal to Iy — I 0 ) at the time 
of switch turn-off gives the magnitude of current which is required to be 
forced turned-off. 

50 
40 
30 
20 

t- 10 

I 0 (Amps.) 

0 

-10 
-20 
-30 
-40 

0 le-05 2e-05 

Time-period (sec.) 

Figure 2.29: Plots of I D to compare with simulation result. 


2.6 Design Equations 

2.6.1 Approximate Design Equations With Sinu- 
soidal Output Current 

The approximate design equations for class-E amplifier can be obtained for 
optimum operating condition. For unity efficiency, input power Pi is equal 
to output power P 0 . Hence, for given power output and supply voltage, 

V 2 

R*=-jT P-58) 

-*< 0 

Where, R dc is equivalent resistance of the amplifier as seen by the power 
supply. Hence, 





(2.59) 
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Using Eq. (2.31) and (2.59), 



(2.60) 


Since, at optimum operating condition, efficiency is unity. Substituting 
this efficiency in Eq. (2.34), we get, 


v, = %? (2.61) 

111 

Hence, substituting V s from Eq. (2.61) in Eq. (2.59) and using Eq. 
(2.29)and (2.30). we get, 

R = 0.5768 R dc (2.62) 

For a given value of operating frequency f, supply voltage V s and power 
output P 0 , the values of load resistance R and shunt capacitance Ci can be 
obtained using Eq. (2.58)-(2.62). 

For a given quality-factor Qi and operating frequency f, output induc- 
tance is given by, 


_ QiR 

27 r/ 

Here, high chosen value of Qi reduces harmonic content in output. 

Using Eq. (2.21), (2.30), (2.60) and (2.62), we can obtain value of C 0 , 
as given below, 



C 0 = 


1 1.1534 

w 2 L 0 1 + <2; — 1.1534 


(2.64) 


The above design equations show that for optimum operation of class-E 
amplifier, there is a unique relation among circuit elements i.e. , there is 
only one set of circuit elements are possible for a given frequency and Qi, 
power output Po and supply voltage Vs. 


2.6.2 Approximate Design Example 

The design example is for an amplifier to deliver power 2 KW at 50 KHz, 
V s = 100V, Li = 5 mH (RF CHOKE), D=|, and a given quality factor Q t . 
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Using equations of previous section, we have obtained approximate val- 
ues of circuit elements. These values are very near to the circuit elements 
which give optimum class-E characteristic. As we increase quality factor, 
the output current becomes more sinusoidal, hence, the fundamental com- 
ponent follow the actual waveform. Thus, the approximate values tend to- 
wards the optimum values with increase in quality factor. Table-2.3, shows 
approximate design values with four value of quality factor. 


Table-2.3 

Approximate Circuit Design Values 

Qi = 5.0 

Qi = 10.0 

Qi = 20.0 

Qi = 40.0 

R = 2.8840ft 

Ci = 0.2026/jF 

C 0 = 0.2869/jF 

L 0 — 45.9003/jH 

R = 2.8840ft 

Ci = 0.2026//F 

C 0 = 0.1247^F 

L 0 = 91.8007/dH 

R = 2.8840ft 

Ci = 0.2026/jF 

C 0 = 0.0586/^F 

L 0 = 183.6014^iH 

R = 2.8840ft 

Ci = 0.2026 /jF 

Co = 0.0284/jF 

L 0 = 367.2028/jH 


Effect of Quality Factor On Circuit Responses 

The four design values shown above in table-2.3 have been simulated. Fig. „ 
2.30 to Fig. 2.32 shows circuit responses V c i, I\ and I 0 for one cycle, with 
different quality-factors. Curves are numbered 1 to 4. Curve 1 is for Qi = 5, 
curve 2 is for Qi = 10, curve 3 is for Qi = 20 and curve 4 is for Qi — 40. 
Fig. 2.30 shows that peak value of V c i goes on decreasing with increasing 
Qi , slope at 6 — 2tt also moves towards zero. 

Fig. 2.31 is zoomed portion of Fig. 2.30 near 6 = 2ir. It is observed 
that with increase in quality-factor, our circuit values tend to be optimum. t 
Fig. 2.32 is plot of I\ and I 0 . It is observed that there is very little decrease 
in Ii with increase in quality-factor. As Ql is increased I 0 at 6 = 2tt tends 
to equal to input current I\. Phase shift in output current is observed with 
quality-factor. 


2.6.3 Modified Procedure for More Accurate Design 

From the harmonic analysis of output current J 0 , it is evident that first 
three harmonics are significant. If in our design equations, three or more 
harmonic components are considered, we will get more accurate values of 
the circuit elements. In this procedure, we use the approximate designed 
circuit element values of previous section as starting point. 






40 


400 
350 
300 
250 
200 

V& (Volts.) 

150 
100 
50 
0 

0 le-05 2e-05 

Time-period (sec.) 

Figure 2.30: Plots of V c i with different Qi. 
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Figure 2.31: Magnified plot of V cU near 9 = 2tt with different Qi- 

We find harmonics of V cl (9), assuming only the fundamental current 
The V c i (9) from Eq. (2.13) can be written in terms of its com- 
ponents, 

V cl (9) = V clav +Vc n (2.65) 

V Op, = 2 [a n cos (n9) + b n sin(n6>)] (2.66) 

n= 1 


Where, k — > co. and 
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Figure 2.32: Plots of I x and I 0 with different Qi. 

1 1 — (— l) n 7 r 

a n = — —[ + -z Iom i cos fa 

TfwCi n z 2 

ll Jomi sin <pi] (2.67) 

1 — nr 

, lrTT 1 ~ (— 1)” 

b n = — —[ — h — —Iorrii cos fa 

ttwCi n n 

-l_ n(l + ( l) _)/ omi cos ^ Iorrii sin 0i] (2.68) 

1 - n 2 2 

At steady-state, the n t/l harmonic output current 7o n , can be found by 
dividing n th harmonic Vcn with n th harmonic series impedance Zo n . Hence, 

Io n 

Where , 

Zo n 
and, 

VC* 
and, 

X n 
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After obtaining, n — harmonic components of Io n , we find new value of 
fundamental current, by putting zero magnitude and zero slope conditions 
of shunt capacitor voltage V cl (6) at 6 = 2i r. V cl {9) is given by Eq. (2.45), 
is used to satisfy the design conditions. Hence, putting V& ( 0 ) equal to zero 
from Eq. (2.45), we obtain; 


1 k 1 — (— l) n 

Iomi cos <f>\ = - - [i rfy + Y" — — I om n cos 6 n ) (2.70) 

2 u 

n—l 


Imposing another design condition — = 0.0 for optimum operating 
condition, we obtain; 


k 

Iomi sin = h — Iom n sin <p n (2-71) 

n= 2 

With modified value of fundamental current, we obtain from Eq. (2.47); 


_ 1 r^ 2 r v—' r r 1 — ( — !) n • , (-l) n ^ , n , n ^ 

Ci = 2™vh h + £ -p— sm cos *->] (2 - 72) 

Since, the value of C\ has been modified, hence, the voltage V c \(6 ) will 
also be modified. New, n th harmonic components of V c \ (6), will depend on 
all harmonic components of output current (not only the fundamental !). 
Hence, Eq. (2.66) will be modified. Suppose modified coefficients are a m 
and b m for m th harmonic of V c i(6), hence, 


OjtTI 


1 r {l-(-l) m } r 7T _ 

r —h + T I°Trim COS (pra 


tvwC\ L m 2 x ' 2m 

+ W — — k 5 — -Iom n sm <fi n \ 


71 = 1 , 71^771 


(2.73) 


and, 


1 ?r 7r A m{(-l) m+n - 1} r , 

1 r 71 T r - i . IV J -/oTO n COS<fy 


-{--h 


n l ix — sin 4> m + r)(m 2 — n 2 ) 

TTwC\ m 2m n=Wm n y n n ) 

| A (-!)"{! -(-l) m > 

71=1 


Iom n cos 4>n] 


mn 


(2.74) 
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From above equations; we find, 


Vc n = y/a^ + arctan (2.75) 

Vn 

Hence, new value of fundamental series impedance will be, 


Zo\ = 


/ of + 6? 

Iomi 


Here onwards, we have two design options: 


1. Design with VI., P 0 , Q t and f constant. 

2. Design with R, Qi and f constant. 


(2.76) 


Design with V s , P 0 , Qi and f constant 


Since, the output current has been modified from Eq. (2.70) and (2.71); to 
keep output power P 0 constant, the new value of output resistance R will 
be, 


R = 


2P 0 


KLi Iomi 


Hence, with fixed Qi series inductor L 0 is given as, 


(2.77) 


L o = 0^ (2.78) 

w 

Hence, new value of series capacitor C 0 is given by, 

Co = ~7 (2.79) 

w 2 Lo - vJ-jZoi - R 2 

a Here, we have new values of R, L 0 , Ci and C 0 ; now we can find new 
values of Zo n . With new values of Zo n and Vcn, we again find Io n from 
Eq. (2.69), and repeat the whole procedure from Eq. (2.69) to (2.79), till 
the values of C ly R, L 0 and C 0 are stablized (i.e., the values computed in 
two consequtive iteration are same). 

Using this procedure, we have computed four designed values, shown 
in table-2.4. The values are with considering only fundamental, only first 
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two, only first three and only first ten harmonics of output current. It is 
evident that as we consider more harmonics, the designed values become 
more accurate, but consideration above three harmonics does not contribute 
significant change in the design values. 


Table-2.4 

Circuit Design Values With More Accurate Procedure 
Design With V s , P 0 , Qi and f constant 


Considering N Number of Harmonics 


N=1 

N=2 

N=3 

O 
! 1 

II 

R = 2.8840C 

Ci = 0.2026/uF 

C 0 = 0.1247/iF 

L 0 = 91.8007/zH 

R = 2.7708Q 

Ci = 0.2275//F 

C 0 = 0.1305/iF 

L 0 = 88.1976/iH 

R = 2.7547Q 

C x = 0.2280/uF 

C 0 = 0.1314^F 

L 0 = 87.6858/jE 

R = 2.74840 

Ci = 0.2292 ^F 

C 0 = 0.1317/iF 

L 0 = 87.4843/iH 


Design with R, Qi and f constant 

In this case, we have fixed R, Qi (i.e., also L 0 ) and f. In the design procedure, 
the output current Io n is changing to satisfy the design conditions of Eq. 
(2.70) and (2.71). Hence, with constant R, the output power P 0 will not be 
constant. The R is hold fixed, at the value designed by approxmate design 
method. At optimum operating condition, the input power should be equal 
to the output power, because of unity efficiency. To satisfy this condition, 
•we will have to vary the input supply voltage V s . m (s 

After Eq. (2.76), we apply following procedure. The new output power 
is given by, 


P n = 


Z k m=l IomlR 


(2.80) 


New supply voltage V s to make input power equal to output power, 


V. 


Po 

h 


(2.81) 


Also, new value of series capacitor C 0 is given by, 


C„ = 


w 2 L 0 - w\JZo \ - R 2 


(2.82) 


Here, we obtained new values of P 0 , V s , C 0 and Vc^. Now, we calculate 
new values of Zo n and find Io n from Eq. (2.69). We repeat the procedure 


fzLyA , 
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from Eq. (2.69) to (2.76) and (2.80) to (2.82), till the values of C x and C 0 
are stablized. 

Using this procedure, we have computed four designed values shown in 
table-2.5, with N=l, N=2, N=3 and N=10. It is evident that only C\ and 
C 0 has been modified. The modified values for N > 3 makes the circuit to 
operate in optimum operating condition. 


Table-2.5 

Circuit Design Values With More Accurate Procedure 
Design With R, Qi and f constant 


Considering N Number of Harmonics 


N=1 

N=2 

N=3 

N= 10 

R = 2.884012 

Ci = 0.2026 /jF 

C 0 = 0.1247^F 

L 0 = 91.8007/xH 

R = 2.8840ft 

Ci = 0.2186/iF 

C 0 = 0.1254/rF 

L 0 = 91.8007//H 

R = 2.8340ft 

Ci = 0.2178/uF 

C 0 = 0.1255/iF 

L 0 = 91.8007/rH 

R = 2.8840ft 

Ci = 0.2184^F 

C 0 = 0.1256/iF 

L 0 = 91.8007/xH 


The flow-chart and computer program of the modified procedure with 
more accurate design is given in appendix (C). 


Circuit Responses With More Accurate Design Values 

Fig. 2.33 to 2.36 shows plots of circuit responses V cl , I x and I 0 with four 
designed circuit values, as it is shown in table-2.4. From Fig. 2.33 curve no.l 
is with design considering fundamental harmonic only, curve no.2 is with 
design considering fist two harmonics, curve no. 3 is with design considering 
first three harmonics and curve no. 4 is with design considering first ten 
harmonics. From the figures, it is evident that only first two harmonics 
are significant. Only curve no.l is distinct, rest of the curves are almost 
covered by curve no. 10. As we include more and more harmonic components 
the design becomes more and more accurate. Fig. 2.34 shows the zoomed 
portion of Fig. 2.33 near 9 = 2tt. In Fig. 2.34 curves are distinct, and it 
can be observed that as we include more number of harmonics, we reach the 
zero magnitude and zero slope condition of V&. Fig. 2.35 shows that output 
current I Q at 9 = 2tt tends to equalize the input current h in magnitude, 
as we consider more and more number of harmonics. Fig. 2.36 is zoomed 
plot of curve shown in Fig. 2.35 near 9 = 0.0. It is clearly observed the 
movement of output current I 0 towards input current Ij at 9 = 0.0. There 
is very small decrease in input current I\ as we consider more number of 
harmonics. This is because, since the output power is constant, the input 
power will be higher than the output power to meet the losses at turn-on. 
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Figure 2.33: Plots of V& with different harmonics in design (P 0 constant). 
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Figure 2.34: Magnified plots of V cl , near 6 = 2% with different harmonics 
in design ( P 0 constant). 
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Figure 2.35: Plot of I x and I 0 with different harmonics in design (P 0 
stant) . 
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Figure 2.36: Magnified plots of I x and I a near 6 = 0.0 with different har 
monies in design ( P 0 constant). 
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Figure 2.37: Optimum circuit response with design assuming R, Qi and f 
constant. 

Fig. 2.37 is the plot of V c i, I\ and I 0 with circuit values given in table- 
2.5 for N=10. The circuit responses show that the operating condition is 
optimum. 




Chapter 3 

Analysis and Design of 
Class-E Push-Pull 
Configuration 


3.1 Introduction 


In the basic single ended class-E amplifier analysed in chapter -2, we have 
observed that the load current and hence load voltage for periods zero to 
DT and DT to T were not symmetrical, even in steady state. This is 
because, during these two intervals the circuit configurations were different. 
Therefore, this single ended class-E may contain even harmonic components 
in the output. For the application as an RF amplifier a second harmonic 
filter is required for this configuration. 

Class-E push-pull amplifier provides a symmetrical configuration as shown 
in Fig. 3.1. It consists of two switches Si and S2 and two shunt capaci- 
tors Ci and € 2 - The Ci and C 2 may include any capacitance inherent in 
switch SI, diode D1 and switch S2, diode D2 respectively. Series tuned 
network consists of L 0 , C 0 and load resistance R. L x and L 2 are high value 
inductances (RF Chokes) in series with voltage sources to realize current 
sources. 

A typical waveforms of push-pull class-E amplifier is shown in Fig. 3.2. 
The voltage magnitude and its slope is always zero across the switches at 
the time of turn-on. At optimum operation, the voltage across Si and S2 
are V cX and V c2 respectively, in switch-off duration. The output current I 0 
is almost sinusoidal. 

The class-E push-pull amplifier can be operated with two or three states 
as shown in Fig. 3.3 and Fig. 3.4. With two modes of operation the duty 
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cycle is 50-percent. For three modes of operation the duty cycle is taken 
greater than 50-percent. 



Figure 3.1: The push-pull class-E amplifier. 



Figure 3.2: Typical sketch of push-pull class-E waveforms Vcl, Vc2 and Io 
at steady-state. 



Figure 3.3: Two modes of push-pull class-E amplifier operation. 
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Figure 3.4: Three modes of push-pull class-E amplifier operation- 
al Exact Analysis of Push-Pull Class-E Am- 
plifier With Two Modes of Operation 
(Using State-Equations) 

The class-E push-pull circuit to be analyzed is shown in Fig. 3.1. The 
present circuit have six energy storage elements, hence, there are six vari- 
ables. The variables Ii, I 2 , I 0 , V&, V c2 and V co are associated with Lx, L 2 , 
L 0 , Cx, C 2 and C 0 respectively. To maintain symmetry in the circuit , the 
values of Cx and C 2 , L\ and L 2 should be same. Also ratings of switches 
Si and S2, diodes Dl and D2 should be equal. Circuit equations have been 
written for two states of operation. 


Circuit Equations for Mode-1 

The equivalent circuit diagram for the circuit mode-1 Si “off’ and S2 “on” 
is shown in Fig. 3.5. Since S2 “on” is associated with diode in antiparallel 
(which conducts for reverse current), it is assumed as a bidirectional switch 
during “on” period. The V C 2 is zero for this interval. There are five first 
order differential equations for this state. 



Figure 3.5: Equivalent circuit for mode-1. 


dli _ V S - Vd 

dt Li 
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dVci _ h - Ip 
dt C\ 

dh = V L 

dt L2 

dh_ _ V* — Ycq _ R t 
dt ~ L 0 ~ T 0 1 
dVco _ h 
dt C 0 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


Circuit Equations for Mode- 2 

For mode-2 the switch state is Si “on” and S2 “off” as shown in Fig. 3.6. 
The five first order differential equations have been written. For this state 
V cl is zero. 



Figure 3.6: Equivalent circuit for mode-2. 


dj 1= v L 

dt L\ 

dh _ Vs ~ V* 

dt L 2 

dV C 2 _ I 2 Io 

dt C 2 

dip _ Vg + Vco 
dt L 0 

dVco dp 

dt C 0 



(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 


Following optimum circuit design values have been taken. 

LI = L2 = 5.0 mH] C1 = C2 = 0.20264 pF; Co = 0.062374/xE 
Lo = 183.601 y.H\ R = 5.768Q; Vs = 100.0V; / = 50 KHz; D = 0.5 
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The simultaneous differential equations for the two states have been 
solved with above circuit values, and the steady-state values of parameters 
(h,h-Jo,V c i,V c2 ,V co ) have been shown in Fig. 3.7 to Fig. 3.10. The flow- 
chart and computer program for solution of the exact analysis with “two 
modes of operation” are given in Appendix- (D). 

Fig. 3.7 is plot ofV c i and V& for one time-period at steady-state. These 
curves satisfy the optimum operating conditions of the push-pull am p lifi er 
The peak value of V cl and V c2 is 358.25V. 

Fig. 3.8 is curve for V co vs time-period. The first half of time-period is 
for SI “off” and S2 “on” and the second half is for S2 “off” and Si “on”. 
The curve is almost sinusoidal and peak magnitude of voltage is 1904V. 

Fig. 3.9 is curve of I x , I 2 , I s i, I s2 for one time-period, since, source 
inductors L x and L 2 and voltage source value V s are equal at the two ends, 
we have I\ = I 2 = 20.1 A. The current source values are almost const an t 
I si and I s 2 are same in magnitude, but 180-degree out of phase. The peak 
switch current J sl and I s2 is 57.4A. For optimum operating condition, when 
switch is off, the current I sX and I s2 flow into capacitors C x and C 2 , respec- 
tively. At the time of switch “off” the magnitude of switch current I s is 
40.01 A. 


Fig. 3.10 is plot of output current I 0 for one cycle of operation. It is 
almost sinusoidal. The peak magnitude of I 0 is 37. 3A. 


Vcl and V c2 

(Volts.) 



Figure 3.7: Curves of V cX and V c2 for one cycle at steady-state. 
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0 le-05 2e-05 

Time-period (sec.) 

Figure 3.8: Curve of V co for one cycle at steady-state. 



0 le-05 2e-05 

Time-period (sec.) 


Figure 3.9: Curves ofh,I 2 Jsi and I s2 for one cycle at steady-state. 


Realizing Current Source as Voltage Source in Series With a 
Large Inductor 


For analysis of push-pull class-E amplifier, the voltage source in series with 
a large inductor (RF Choke) can be replaced by a current source as shown 
in Fig. 3.11. The current source magnitude h is decided by the magnitudes 
of voltage source and effective impedance in series with it. The power input 
to the amplifier is decided by the magnitudes of I\ and J 2 . The replacement 
of voltage source and large series inductor by a current source, makes the 
analysis simple, because the number of state equations are reduced. 





38 
28.5 
19 
9.5 

(Amps.) 0 

-9.5 
-19 
-28.5 
-38 

0 le-05 2e-05 

Time-period (sec.) 

Figure 3.10: Curve of I 0 for one cycle at steady-state. 




Figure 3.11: Current source realized by voltage source in series with a large 
inductor. 

Voltage source in series with large inductor was replaced by the current 
source in Fig. 3.1, and the circuit was simulated. It gave the same results 
as shown from Fig. 3.7 — 3.10 for the same values of input currents and 
circuit elements. 


3.3 Harmonic Analysis of Output Current 


The push-pull class-E amplifier may contain harmonics in its output current 
wave. Because of symmetry, the even harmonics will be absent. Using 
“Discrete Fourier Transform Method”, the harmonic components of output 
current I 0 shown in Fig. 3.7, have been obtained. The peak magnitudes 
and phase angles are tabulated in Table-3.1. The harmonic components 
were calculated using NAG subroutine C06EAF. 
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Table- 3.1 

Harmonics and Phases of Output Current 

Harmonic 

Peak magnitude (Amp.) 

Phase angle (Degree) 

1 

37.105 

-31.996 

3 

0.2904 

183.261 

5 

0.05544 

181.236 

7 

0.01891 

176.874 

9 

0.008546 

168.236 

11 

0.004814 

154.744 

13 

0.003382 

139.128 


From above table, it is evident that magnitudes of third and above har- 
monic components are negligible in comparison to fundamental component. 
Hence, the output current waveform is almost sinusoidal. 


3.4 Analysis Assuming Constant Input Cur- 
rent and Sinusoidal Output Current With 
Two Modes of Operation 

3.4.1 Analysis 

■For approximate analysis and design, the equivalent circuit of class-E push- 
pull amplifier is shown in Fig. 3.12. To analyze the circuit following as- 
sumptions have been made. 

• Current through the load I 0 (6) is pure sinusoidal. 

• Switches SI, S2 and diodes Dl, D2 are ideal. 

• Capacitors Ci and C 2 may include any capacitances inherent in Si 
and S2 respectively and are independent of any voltage across the 
switches. 

• Inductors L x and L 2 are large enough to make currents h and I 2 
almost constant. 

• Li and L 2 are equal. Switches SI and S2, diodes Dl and D2 are of 
same ratings for symmetrical operation of the circuit. 

For optimum operation of class-E push-pull amplifier the switch oper- 
ating frequency will be different from the resonant frequency of L 0 and C 0 . 





Figure 3.12: Equivalent circuit of push-pull class-E amplifier for approxi- 
mate analysis. 

Another equivalent circuit shows, series tuned load circuit consisting of jX, 
R and a voltage source V c , as shown in Fig.3.13. The reactance X is equal 
to ( wL 0 — ^-). Resistance R includes any series resistance inherent in L 0 , 
C 0 and equivalent ‘lumped’ a.c. resistances of switches and diodes. The 
voltage source V c {6) is equal to V ci (6) from 0 to 7r and is equal to —V& (0) 
from 7 r to 2rc. The waveforms of V e i(6) and V c2 (8) are shown in Fig. 3.2. 
The voltage source V c (9) is alternating, hence, output current I 0 is also 
alternating. 



Figure 3.13: Equivalent circuit of push-pull class-E amplifier for approxi- 
mate analysis. 

For the case of Si “off” and S2 “on”, the equivalent circuit given in Fig. • 
3.13 will be applicable form 0 to 7r, and for case of SI “on” and S2 “off”, 
the equivalent circuit will be applicable from 7r to 2tt. The output current 
and voltage are assumed sinusoidal. Thus, 

I 0 (9) = Iom sin (0 + <f>) (3.11) 

V o (0) = IomRsin(e + <£) (3-12) 

Where, 0(=wt) is “angular time”, peak magnitude I om and phase angle 
is <fi (as shown in Fig. 3.2), are to be determined. V C {Q) is sum of the 
voltages V x {6) and Vr{9). Hence, 
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Vc{0) = Iom^/(R 1 2 + X 2 ) sin(# + 4> + ib) (3.13) 

Where, ip = arctan (=|). 

V c (6 ) is represented by following equation; 

V c (e) = V cl (6) - V c2 (e) ( 3 . 14 ) 

Where, V& (6) and V c2 (6) can be obtained by following method. 

For the duration Si “off” and S2 “on” V c2 (6) is zero, and 

1 r e 

Vci(0) = J o [-^1 - hm sin(u 4- (p)]du (3.15) 

Hence, 

V*(6) = + lom cos (6 + (p) - Iom COS <p] (3.16) 

Where, 6 is between 0 and n for D= 0 . 5 . 

For the duration SI “on” and S2 “off” the V cl (9) is zero, and 

V c2 {6) = £ [I 2 + I om sin(u + 4>)]du (3.17) 

Hence, 

V c 2 (9) = -^r[h(0 -7 r) - I om cos (8 + <p) - I om cos <p] (3.18) 

wU 2 

Where, 6 is from 7 r to for D=0.5. 

The average voltage across inductor L\ or L 2 is zero. Hence the voltage 
across the source is equal to the average voltage across capacitors Ci or C 2 . 
Therefore, source voltage, 

1 r 27r 

V, = — / V cl me (3.19) 

Z7T J 0 

From equations (3.16) and (3.19), we have, 

1 7 r 2 

= n 7T hr A - 2 Iom sin <f> - Ttlom COS <j>] 

2i7rwC\ 2 


(3.20) 



59 


Also, the source voltage, 


v ’"hC Va{9)dB (3 - 21 > 

From equations (3.18) and (3.21), we have, 

1 7T^ 

Vs = 2 ^'g f Y^ 2 ~ 2 ^ om sin ^ cos 0] (3.22) 

From equations (3.20) and (3.22), it is evident that for symmetry, we 
should have equal voltage sources V s , equal shunt capacitors Cj and C 2 , 
equal RF Choke inductors L\ and L 2 . The voltage across the series impedance 
Z 0 , is V c {6). Where, V c (9) is given by Eq. (3.14). Fundamental component 
of V c (6) is obtained as, 


Vc fundamental ( d ) = ^t-[(-4/i+7 r/ om cos <f>) cos 9+ {2ttI v -7rI om sin (j>-AI om cos ©)] sin 6 

(3.23) 

In above equation, J 2 and C 2 have been replaced by I\ and C\ respec- 
tively. From Eq. (3.13), V c (9 ) can be written as, 


V c (9) = [. RIom sin 6 + XI om cos 4>] cos 9 + [RI om cos 4> - XI om sin <fi] sin 6 

(3.24) 

Equating coefficients of sin# and cos 9 from Eq.(3.23) and (3.24), we 
get following equations. 


— AIi ~ nwCiRIom sin (f> + 7r(l — wCiX)I om cos (j) = 0 (3.25) 

and, 

2t rrli + tt(-1 + wCiX)I om sin <\> - (4 + txwC x R )I om cos tj> = 0 (3.26) 

Using Eq. (3.25), (3.26) and (3.20), we get following matrix equation; 


—4 —ttwCiR (7t — nwCiX) 


‘ h 


' 0 

2tt (ttwCiX — 7 r) (—4 — ttwCiR) 


lom Sin 0 


0 

-T "2 


lorn COS (j) 


27r7uCiVs _ 
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From above matrix equation, we can solve for I lt I om and <j>, for a given 
set of circuit element values. The angle <f> is obtained as, 

= arctan [ — t ] (3.27) 

W 2 wCi J 

Analysis For Optimum Operation 

Since, all circuit elements are ideal, there are no losses. The loss occurs 
only if the switch-on take place at finite V c \ or V c2 . The requirement for 
a hundred-percent efficiency can be realized by setting equation (3.16) and 
(3.18) equal to zero at 0 = -it and 0 = 2n respectively. Therefore, 

Vci {9)/e=ir = o = I id + J om cos(# + 4>) - I om coscp (3.28) 

Which gives, 

cos (f> = (3-29) 

The slope of V cl (0) at the time of turn “on” of Si must also be zero for 
optimum operation. Hence, 

dV c i(0) ^ _ o = J x - J om sin(d + 4 >) (3.30) 

at 

Which gives, 

sin^> = — y~ (3.31) 

-Lorn 

From Eq. (3.29) and (3.31), we get, 

<f> = -32.481637° = —0.5669115° (3.32) 

From Eqs. (3.20), (3.29) and (3.31), we get 

Ii = Vg-irwCi (3.33) 

Similarly, because of symmetry of the circuit following can also be ob- 

tained. 
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h = V s 7 twC 2 (3.34) 

From Eq. (3.31) and (3.32). 

Iom = 1.8621/x (3.35) 

Peak switch current is given by, 

Islpk = "t Iom (3.36) 

From Eq. (3.35) and (3.36), we get 

Isipk = 2.8621 I x (3.37) 

For the value of 9 vmax at which the value of V c i (9) is maximum, from 
Eq. (3.16); 

— ~ - / e=e vma x = 0 = Jx — I om sin (9 vmax + 6) (3.38) 

Using Eq. (3.31), (3.32) and (3.38), we get 

Ovmax = 64.9633° = 1.133823 c (3.39) 

Substituting value of 9 vmax in place of 9 in Eq. (3.16) and using Eq. 
(3.31), (3.32) and (3.33), we obtain, 

Kw = 3-562U, (3.40) 

Similarly, we can also obtain, 


Uc2 


max 


3.562V; 


Output power is given by 


P 0 = 


Safi 

2 


(3.41) 


(3.42) 


Input power is, 
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Pi = V s (I 1 + I 2 ) 

Efficiency of the amplifier is. 


ILR 

2V s (h + I 2 ) 


(3.43) 


(3.44) 


3.4.2 Comparison of Approximate Analysis Results 
with Results of Simulations 

The results obtained by above analysis will be approximate because the 
output current is assumed sinusoidal and input current is assumed co ns t an t 
For a given value of V s = 100V, / = 50 KHz and C x = C 2 = 0.20264 /j,F, we 
calculate I\ and I 2 20A from Eq. (3.33) and (3.34). The simulation result 
in Fig. 3.9 shows its value 20. 1A. 

The peak output current I om given by Eq. (3.35) is calculated as 
37.242A. The simulation result in Fig. 3.10 shows 37.3A. 

Peak switch current is calculated as 57.242A from Eq. (3.37) and 
simulation result gives 57.4A, as in Fig. 3.9. 

The values of V c i max and V C 2 max are calculated as 356.2V from Eq. (3.40) 
and (3.41). The simulation result gives V c i max = V c2max = 358.25V from 
Fig. 3.7. 


3.5 Design of Push-Pull Class-E Amplifier 
For Optimum Operation at D=0.5 

3.5.1 Design Equations 

For symmetrical operation of circuit, we have to keep C\ = C 2 . The supply 
■voltage is common to both ends. Also, Li and L 2 are kept equal. From 
Eq. (3.33) and (3.34), it is observed that Rdcl(= ^ and Rdc2(= 

V- = — V-) are the resitive loads that the amplifier shows to respective 
power supplies at the two ends(Fig. 3.1). For Rdcl=Rdc2=Rdc and unity 
efficiency, we have, 


Rdc = 


2 Vi 

Pc 


(3.45) 
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Where, P 0 is equal to sum of input powers at both end of the circuit at 
optimum operating condition. From Eq. (3.44) with I x = I 2 and T] = 1 , we 
get, 


T 2 r> 
t/ 1 om rL 

4/x 


(3.46) 


Hence, substituting from Eq. (3.33) and using Eq. (3.31) and (3.32) 
(also Rdc = we have, 

R = 1.1536i?dc (3.47) 

Using Eq. (3.47) and Rdc = we have, 


„ _ 0.0584421 

Cl ~ Jr 

For a given quality factor Qi, 

_ QiR 
°~ 2t rf 


(3.48) 


(3.49) 


Substituting calculated values of R, C u f, L 0 (from above equations) 
and (f> (from Eq. (3.32)) in Eq. (3.27), we obtain following equation for C 0 \ 


1 1.1524959 

w 2 W l + Q t - 1.1524959 


(3.50) 


Using Eq. (3.45) — (3.50) the circuit elements can be designed for op- 
timum operating condition. The value of L\ and L 2 should be chosen as 
large as possible to maintain input current almost constant. 


3.5.2 Design Example 

Same design problem has been taken as it was in chapter-2, with double 
amount of power. The amplifier is to deliver 4KW at 50 KHz,D=0.5, V s = 
100U, L x = L 2 = 5.0mH (RF Choke), and Q t = 10. 

Using design equations of previous section the designed values obtained 
are given below: 

C X = C 2 = 0.20264 nF; C Q = 0.062374 /nF; L 0 = 183.601^77; R = 5.7680 


The push-pull class-E amplifier operates in optimum operating condition 
with the above designed circuit element values. 
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3.6 Exact Analysis of Push-Pull Class-E Am- 
plifier With Three Modes of Operation 
(Using State-Equations) 

The push-pull class-E amplifier shown in Fig. 3.1 is to be analyzed with 
three modes of operation (shown in Fig. 3.4). In this case also, the 
symmetry of the circuit has been maintained. There are six variables 
-f i : -f 2 j Jo, V c i, V C 2 and V co associated with six energy storage circuit elements 
L\, L 2 , L 0 , Ci, C 2 and C 0 respectively. Circuit equations have been written 
for three modes of operation. 


Circuit Equations for Mode-1 

Equivalent circuit diagram for mode-1 is shown in Fig. 3.5. For this mode 
switch Si is “off” and S2 is “on”. During “on” period the switch S2 conducts 
forward current and diode D2 conducts reverse current. The voltage V C 2 is 
zero for this mode. There are five first order differential equations for this 
state. 


dh _V S - V cl 
dt Li 

(3.51) 

dV cl h - 1 0 
dt Ci 

(3.52) 

dh_V L 
dt L 2 

(3.53) 

Ki - F co R 

L 0 L 0 0 

(3.54) 

dVco do 
dt C 0 

(3.55) 


Circuit Equations for Mode-2 

For mode-2 the switches Si and S2 are in “on” state as shown in Fig. 3.14. 
Due to presence of antiparallel diodes, the switches will be conducting in 
bidirection. For this case both V& and V C 2 are zero. There are four first 
order differential equations for this mode. 
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Figure 3.14: Equivalent circuit for mode-2 in three mode operation. 


dh _V L 
dt Li 

(3.56) 

dh_V L 
dt L 2 

(3.57) 

dl 0 V co R 

dt ~ L 0 L 0 ° 

(3.58) 

-?lc? 

II 

kS -+0 
"C3 

(3.59) 

Circuit Equations for Mode-3 


For mode-3 the switch state is Si “on” and S2 “off”, it is same as mode-2 of 
“two mode operation” as shown in Fig. 3.6. The five first order differential 
equations have been written. For this state V cl is zero. 

dh_Vs_ 
dt Li 

(3.60) 

<3 

C4 

i ^ 

II 

1 ^ 

(3.61) 

dV c 2 I*l + Io 
dt - C 2 

(3.62) 

14 

i 

E? 

1 ° 
+ ^ 

£ 

II 

(3.63) 


dVgo I O 

dt C 0 


(3.64) 
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Following optimum circuit, design values have been taken. 


LI = L2 = 5.0771#; Cl = C 2 = 0.03168/uF; Co = 0.02839 fxF 
Lo = 381.0682/nF/; R = 11.971GH; Vs = 100. OF; / = 50 KHz; D = § 

For three modes of operation the duty cycle will be greater than 50-percent. 
For present case the duty-cycle is chosen as §. The circuit operates for 
| of time-period in mode-1, then ± of time-period in mode-2, then | of 
time-period in mode-3 and again in mode-2 for | of time-period, then 
a, gain in mode-1 for | of time-period, like this the cycle repeats. The 
simultaneous differential equations lor the three-modes have-been solved 
with above circuit values, and the steady-state values of circuit parame- 
ters (I\, 1‘2, 1»,V c \,Vr2,Vc„) are shown in Fig. 3.15 to 3.18. The flow-chart 
and computer program for solution of the exact analysis of three modes of 
operation are given in Appendix-(E). 

Fig. 3.15 is plot of. V c \ and V c2 lor one time-period at steady-state. The 
peak values of V c \ and V& is 549.12V. Switch turn-on is achieved at zero 
magnitude and zero slope of F c] and V c2 . Fig. 3.16 is curve of voltage 




Figure 3.15: Plots of V ci and V c2 in three mode operation. 

V co for one time-period at steady-state. This is the voltage across series 
capacitor C () . It is almost sinusoidal with peak magnitude of 2862.2V. 

Fig. 3.17 is plots of currents I\, I 2 , /si, and I s2 for one time-period at 
steady-state. Curents L and I 2 are constant at 20.2A. Current I s y starts 
with 15.87A, this is its I slof/ value. At switch turn-on time, value of I sl 
is zero. The peak value I s i pea k is 46.34A. The current I s2 is same as I sl in 
magnitude but opposite in phase. 
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Fig. 3.18 is plot of output current I 0 for one cycle at steady-state. It is 
almost sinusoidal with its peak magnitude of 26.22A. It starts with a value 
4.36A with positive phase angle. 
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Figure 3.16: Plot of V co in three mode operation. 
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Figure 3.17: Plots of 7i, J 2 ,/sl and I s2 in three mode operation. 




3.7 Harmonic Analysis of Output Current 
For Three-Modes of Operation 


The push-pull class-E amplifier with three modes of operation retains sym- 
metry in its output current waveform. Hence, because of symmetry, the 
even harmonics will be absent. Using “Discrete Fourier Transform Method” 
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Figure 3.18: Plot of I 0 in three mode operation. 



(NAG subroutine C06EAF), the harmonic components of the output cur- 
rent waveform I 0 (shown in Fig. 3. 18) are obtained. The magnitudes and 
phase angles are tabulated in Table 3.2. 


Table-3.2 

Harmonics and Phases of Output Current 

Harmonic 

Peak magnitude (Amp.) 

Phase angle (Degree) 

1 

25.839 

10.77 

3 

0.5497 

226.26 

5 

0.05708 

194.50 

7 

0.02819 

189.89 

9 

0.01059 

199.80 

11 

0.00506 

184.32 

13 

0.00341 

189.48 


From Table-3.2, it is evident that third and higher harmonics are very- 
small in magnitude in comparison to fundamental harmonic component. 
This indicates that the output current I Q is almost sinusoidal. 
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3.8 Analysis Assuming Constant Input Cur- 
rent and Sinusoidal Output Current With 
Three Modes of Operation 

3.8.1 Analysis 

For analysis of push-pull class-E amplifier (shown in Fig. 3.1) with three 
modes of operation, all assumptions are same as mentioned in section 3.4.1. 
Since, the operating frequency is different from the resonant frequency of L 0 
and C Q , the equivalent circuit shown in Fig. 3.12 is valid for this analysis. 
The output current I 0 is assumed sinusoidal. Consider, 


I 0 {9) = I om sin(0 + <f>) 


(3.65) 


Output voltage is, 


Vr{ 6) = IoraR sm(9 + <f>) (3.66) 

Here, 0(=wt) is “angular time” and q> is phase angle. For three modes 
of circuit operation, the equivalent circuit shown in Fig. 3.13 is assumed. 
The series tuned circuit consisting of R, L 0 and C 0 is excited by a voltage 
source value V c (9). Where, V c {9) for different modes (Fig. 3.4) is defined 
as following: 


V c {9) = V c \ (6) from 0 to 9\ for mode-1 
= 0 from 9i to 7r for mode-2 
= -V C 2 (9) from n to (re + $i) for mode-3 
=0 from (7r + 6h) to 2 n for mode-2 

Hence, 

V c (9) = V cl (9)-V d2 (9) (3.67) 

V c (9 ) can also be expressed by; 


V c {9) = 7 om i/(R 2 + X 2 ) sin(0 + f> + VO 


(3.68) 
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Where, ip = arctan (^). 

For the duration 0 to 8 X Si “off” and S2 “on” (mode-1) V c 2 {d) is zero, 

and 

1 r e 

Vci(0) = -^=r J ^ [I i - I 0TO sin(u + 6 )]du (3.69) 

Solving above, 

V cl (#) = + 4m cos (^ + <t>) - lorn cos 4>] (3.70) 

Where, 9 is between 0 and 6 \. 

For the duration k to (^ + 6*].) Si “on” and S2 “off’ (mode-3) the V cX {9) 
is zero, and 

Va (9) = J [I 2 + / om sin( , u + q>)]du (3-71) 

Solving above, 

V c 2 (d ) = -±r[I 2 (0 - 7T) - lorn COS (8 + <!>)- Iom COS <f>] (3.72) 

WO 2 

Where, 8 is from 7r to (7r + 61 ). 

The average voltage across shunt capacitors C\ or C 2 is equal to source 
voltage Vs, because average voltage across L x or L 2 is zero. Therefore, 
source voltage, 

K = i f'v a (S)dB (3.73) 

2n Jo 

From equations (3.70) and (3.73), we get, 

V s = — + I ora sin(^i + — Iom SHI 4> @ l lorn C °S (f)] (3.74) 

2'kwCi 2 

Similarly, the average voltage across C 2 is given by, 

1 9 2 

= — bf 4 + 4m sm(0! + $ - 4m sin <p - 6 X 4m COS <f>] (3.75) 

2'kwCi 2 
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For symmetry, we assume, in the following, equal voltage sources V s . 
equal shunt capacitors and C 2 , equal RF Choke inductors h and L 2 at 
the two ends of the circuit. The voltage across the series impedance Z 0 . 
is V c {9). Where, V c (9) is given by Eq. (3.67). Fundamental component of 
V c (9) is obtained as, 


Vc fundamental ( 9 ) = cos 9 + b x sin 9 (3.76) 

Where, 

1 r 8 1 r(ir+ 6 1 ) 

ai = -[ / y cl (9) cos 9d9 - / V c2 (9) cos 9d6] (3.77) 

7f JO Jir 

and, 

1 r& i /•(w+Si) 

h = -[ V el ( 9 ) sin - / 1^(0) sin 0d0] (3.78) 

Tt JO Jir 

Where, V c i (9) and V c2 (9) are substituted from Eq. (3.70) and (3.72), 
respectively. Solving Eq. (3.77) and (3.78), we get, 


ai = — ^—-\2Ii(dism9i + cos9i-l) + (9i + ]-sm29i-2smdi)I om cos(f> 

ttwCi 2 

+| (cos 2<9 X - 1)1 cm sin <j>] (3.79) 

and, 

bi = — [2Ii(— 9i cos 0i + sin$x) + (2 cos#! — - cos2$x — -)J om cos^> 
ttwCi 2 z 

+(i sin 29 x — 9i)I om sin 4>] (3.80) 

Simplifying the Eq. (3.68), we get, 


VJ9) = (RI om sin <f> + XI om cos <f>) cos 9 + {RI om cos <p - XIom sin <t>) sin 9 

(3.81) 

Equating coefficients of cos 9 and sin$, from Eq. (3.76) and (3.80), and 
substituting values a x and b\ from Eq. (3.79) and (3.80), we obtain following 
equations: 
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2/i ( 6 X sin 6 X + cos 6 X - 1) + (8 X + I sin 26 x - 2 sin 6 X - nwCiX)^ cos <j> 
+{|(cos20i ~ 1) - nwClR} I om sin <p = 0 (3.82) 


and, 


2Ii(—8i cos 6 X + sin 9 X ) + (2 cos 0\ — \ cos 20 x — | — irwC x R)I 0Tn cos cj> 

+(7rwC x X + \ sin 20i - 9 x )I om sin <j> = 0 (3.83) 


Using Eq. (3.74), (3.82) and (3.83), we form the following matrix equa- 
tion: 


'T n X12 Xi 3 ' 


‘ h 


' 0 

X21 X22 X23 


Iom sm (j) 

= 

0 

. Xz x X 32 X33 _ 


Iom COS (j) 


2irwC x V s 


Where, 


Xu = 2 (6 X sin 6 X + cos 9 X — 1); X12 = |(cos 26 x — 1) — mvClR 
X13 = 0i + | sin 20i — 2 sin 0i — ttwC x X 
X21 = 2(— 0i cos 0i + sin 0i); X 2 2 = 7rwC x X + \ sin 20i — 0i 
X23 = 2 cos 0i — 5 cos 20i — | — ttwC x R 

X31 = X32 = cos 0 i — 1 ; X33 = sin 0 i — 0 i 


From above matrix Eq. (3.84), we can solve for I x , I om and <j>, for a given 
set of circuit element values. The angle <j> is obtained as (with 9 X = ■y), 


<j> — arctan 


1.8191374 - 19.154173wCiX + nwCRR 
2.0453143 + 19.154173u>(7i-R + %wC x X 


(3.85) 


Analysis For Optimum Operation 

For opt im um operation of the circuit the magnitudes and slopes of V cX (0) 
and U c2 (0) are zero at switch turn-on, i.e., at 0 = 0i and 0 = (n + 0i), 
respectively. For D — |, the angular angle 0 X = Therefore, 

Vd ( 0 ) /e^-f = 0 = h& + hm COS (0 + Iom cos <f> (3.86) 
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Which gives, 

2n 

—Ii - 1.5 lom cos cf> - 0.8660254/ om sin <f> = 0 (3.87) 

The slope of V c \{0) at the time of turn “on” of Si must also be zero for 
optimum operation. Hence, 

dVAB) , 

^ / 0 =?z — 0 — Ii — I om sin(f? + <fi) (3.88) 

Above equation can be written as, 


h ~ 0.8660254/ om cos <f> + 0.5 I om sin <$> = 0 (3.89) 

Eq. (3.87) and (3.89) give following equations: 


cos(j> = 1.275482- y- 

J-om 

and, 


sirup = 0.2091995 


h_ 

1-om 


(3.90) 


(3.91) 


The above two equations yield, 


4> = 9.3144939° = 0.1625685 c (3.92) 

Because of symmetry, we get the same phase angle <f> by applying zero 
magnitude and zero slope condition on V&tf). 

From Eqs. (3.74), (3.90) and (3.91), we get 


V, _ 0-3126826 . qon 

h 2 -kwCi ( ‘ } 

Similarly, because of symmetry of the circuit, following can also be ob- 
tained. 


V, _ 0.3126826 
I 2 2ttwC2 


(3.94) 
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From Eq. (3.91) and (3.92). 


Iom = 1.292524/i 
Peak switch current is given by, 


(3.95) 


^slpk -^1 4” ^om (3.96) 

From Eq. (3.95) and (3.96), we get 


I slpk = 2.292524/x (3.97) 

For the value of 8 vmax at which the value of V c i (9) is maximum, from 
Eq. (3.70); 


^ -/ e=e vmax — 0 — Ii — I om sin (9 vmax + <f>) (3.98) 

Using Eq. (3.91), (3.92) and (3.98), we get 

9 vmax = 41.371° = 0.72206 c (3.99) 

Substituting value of 9 vmax in place of 9 in Eq. (3.70) and using Eq. 
(3.90), (3.91) and (3.93), we obtain, 


Kw = 5. 33517, (3.100) 

Similarly, we can also obtain, 

V a™ = 5.335^, (3.101) 


Output power is given by 


Po = 


T 2 R 

• L om lb 


(3.102) 


Input power is, 


Pi = V,(h + h) 


(3.103) 
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Efficiency of the amplifier is, 


V = 



2V,{I 1 + I 2 ) 


(3.104) 


3.8.2 Comparison of Approximate Analysis Results 
with Results of Simulations 


The results obtained by above analysis will be approximate because the 
output current is assumed sinusoidal and input current is assumed constant. 
For a given value of V s = 100V, / = 50 KHz and C x = C 2 = 0.03168 fj,F, we 
calculate I x and I 2 20A from Eq. (3.93) and (3.94). The simulation result 
in Fig. 3.17 shows its value 20. 2A. 

The peak output current I om given by Eq. (3.95) is calculated as 25.85A. 
The simulation result in Fig. 3.18 shows 26.22A. 

Peak switch current I s Xp k is calculated as 45.85A from Eq. (3.97) and 
simulation result gives 46.336A, as in Fig. 3.17. 

The values of V cXmax and V c2max are calculated as 533.5V from Eq. 
(3.100) and (3.101). The simulation result gives V c i max = V c2max = 549.12V 
from Fig. 3.15. 


3.9 Design of Push-Pull Class-E Amplifier 

cy 

For Optimum Operation at D — 3 


3.9.1 Design Equations 


For symmetrical operation of circuit, we have to keep C\ — (A and L\ — L 2 . 
The supply voltage is common to both ends. Net resistive load offered by 
the amplifier at each supply end is, 



(3.105) 


Where, P 0 is desired output power at unity efficiency. 


Also, 


Rdc 


V 

h 


(3.106) 
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From Eq. (3.104) with I\ — Io and rj = 1, we get, 





(3.107) 


Hence, substituting V s from Eq. (3.106) and using Eq. (3.91) and (3.92), 
we have, 

R = 2.39432267?^ (3.108) 

Using Eq. (3.93) and (3.106), we have, 


c _ 0.3126826 
1 27:wR dc 

For a given quality factor Qi, 

j _ QiR 

L / 0 — 

W 


(3.109) 


(3.110) 


Substituting calculated values of R, C 1; f, L 0 (from above equations), 6 X = 
Y and <j> (from Eq. (3.92)) in Eq. (3.85), we can obtain following equation 
for C Q \ 


1 . 0.6330537 

v?L 0 11 + Q, - 0.6330537 


(3.111) 


Using Eq. (3.105) — (3.111) the circuit elements can be designed for 
optimum operating condition. 


3.9.2 Design Example 

The design problem of Sec. 3.5.2 has been repeated here. The amplifier 
is to deliver 4KW at 50 KHz, 7} = §, = 100V, L x = L 2 = 5.0 mH (RF 
Choke), and Qi = 10. 

Using design equations of previous section the designed values obtained 
are given below: 


Ci = C 2 = 0.03168^F; C 0 = 0.0283 9fxF 
L 0 = 381.0682 nH\ R = 11.97160 

The class-E push-pull amplifier operates in optimum operating condition 
with the above designed circuit element values for three modes of operation. 
This design has been simulated and the results have been shown in Sec. 2.6. 



Chapter 4 


Conclusion 


Based on design and analysis carried out, following conclusions can be 
drawn: 


• Simulation study of basic class-E amplifier shows that, during tran- 
sient period the shunt capacitor voltage V cl and out put current I 0 
increases exponentially and reaches the steady-state in around 400 
cycles. At steady-state, peak voltage across the switch is 3.587 times 
the supply voltage. Peak switch current is 2.8 times the input cur- 
rent from the source. For optimum operation, at the time of switch 
turn-on both the voltage across switch and the slope of the voltage is 
zero. At the time of switch turn-off the switch current is around 2.16 
times the source input current. 

• Power output of the class-E amplifier follows the square law with 
voltage source, this is because, an equivalent resistance Rd c can re- 
place the amplifier across the supply voltage at optimum operating 
condition. 

• If the load resistance is decreased from its optimum value, the power 
output to the load is decreased, and the circuit operates in under- 
damped suboptimum operating condition. Increase in the resistance 
makes the circuit operation in overdamped suboptimum operating 
and power output decreases. The power output is maximum only in 
optimum operating condition. 

• The basic class-E amplifier, when it is used for high frequency induc- 
tion heating, generally the load resistance increases and load induc- 
tance decreases with temperature. For a given characteristic of change 
of R and L 0 with temperature, the circuit can be designed to oper- 
ate in underdamped suboptimum operating condition. Underdamped 
suboptimum operating condition also has 100% efficiency. 
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• For underdamped suboptimum operating condition, the voltage across 
the switch will be much larger than that in optimum operating con- 
dition. 

• Since class-E amplifier consists of tuned R- L 0 - C 0 circuit, a small 
change in operating frequency causes large change in output power. 

• Harmonic analysis of the basic class-E amplifier load current and volt- 
age across the switch, shows that all harmonics are present. At opti- 
mum operating condition the second and third harmonics of output 
current are 5.8% and 0.86% of the fundamental. Only second har- 
monic is significant. 

• An analysis assuming constant input current and sinusoidal output 
current has been carried out for solutions of basic class-E amplifier 
for a given set of circuit elements and for any duty-cycle. 

• Analysis assuming constant input current and sinusoidal output cur- 
rent for optimum operating condition at 50% duty-cycle has been 
done. The analysis assuming constant input current and sinusoidal 
output current is an approximate analysis, but this gives design equa- 
tions to obtain approximate values of circuit elements. 

• Analysis assuming constant input current and non-sinusoidal output 
current has also been carried out. This analysis includes finite num- 
ber of harmonic components in the output current, leading to more 
accurate analysis. For a given set of circuit elements the accurate 
solutions can be obtained for any duty-cycle. 

• A comparison of the analysis assuming sinusoidal output current and 
the analysis assuming non-sinusoidal output current with the exact 
analysis has been done. The analysis assuming non-sinusoidal output 
current gives results very close to those obtained by the exact analysis, 
using state equations. 

• For a given operating frequency, output power, supply voltage and 
quality-factor the approximate circuit element values can be obtained 
for D=0.5. A comparative study of different quality-factor on cir- 
cuit responses shows that, if we choose higher quality-factor, the ap- 
proximate designed circuit values become more accurate. This is be- 
cause, with higher quality-factor the output current becomes more 
sinusoidal. 

• A modified procedure for more accurate design for any quality-factor 
has been contributed. This has two design options: 

1. Design for a given P 0 , Qi and f 

2. Design for a given R, Qi and f 
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• “Modified procedure for more accurate design” assumes finite num- 
ber of harmonic components in the output current. Considering dif- 
ferent number of harmonics (in output current), the circuit design 
values “With More Accurate Procedure” were obtained. A compar- 
ative study of the circuit responses with the different circuit values 
have been done. It shows that, assuming first and second harmonics 
in output current is sufficient to obtain accurate optimum circuit el- 
ement values, but as more number of harmonics were considered in 
the output current, more accurate design values were obtained. 

• Class-E push-pull amplifier have analysis and design procedure similar 
to that of the basic class-E amplifier. The push-pull class-E amplifier 
operating in optimum condition with 50% duty-cycle has peak voltage 
across the switch 3.58 times the supply voltage. Peak switch current 
is 2.86 times the input source current. The current through the switch 
at the instant of switch-off is twice the input source current. 

• Harmonic analysis shows that all even harmonic components are ab- 
sent in the output current. Third and higher order of harmonic com- 
ponents are negligible in comparison to fundamental. Third harmonic 
component is 0.78% of the fundamental. 

• Analysis assuming constant input current and sinusoidal output cur- 
rent with D=0.5 has been carried out. Results of this analysis and 
results of the exact analysis using state equations were compared, 
and both results have almost the same value. The analysis assum- 
ing constant input current and sinusoidal output current for optimum 
operating condition gave equations for the design of circuit element 
values. 

• An analysis and design of push-pull class-E amplifier was carried out 
with D = |. This has three modes of operation. At optimum operat- 
ing condition peak voltage across the switch is 5.49 times the supply 
voltage and peak switch current is 2.29 times the input source cur- 
rent. Power capability of the switch have been reduced, but the switch 
current at the instant of turn-off is only 0.79 times the input source 
current. 

• Harmonic analysis of load current shows that even harmonics are ab- 
sent and third harmonic is only 2.1% of the fundamental. 

• Analysis assuming constant input current and sinusoidal output cur- 
rent with D > 0.5 have been done. This analysis for optimum oper- 
ating condition (with D — §) leads an accurate design of the circuit. 
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Suggestions For The Future Work 

1. A study can be carried out for implementing the class-E amplifier 
as DC to DC converter by changing the duty-cycle (at an operating 
frequency) such that the circuit operation remains either in optim um 
or in underdamped suboptimum. 

2. A study can be done to bring the class-E amplifier operation back to 
optimum by changing either duty-cycle or operating frequency, if due 
to small change in circuit element values, the optimum operation is 
disturbed. 
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Appendix A 


Flow-Chart For Exact Analysis 
of Class-E Amplifier (With 
Diode), Using State Equations 


( START ) 

A 

READ 
Supply Voltage = Vs 
Operating Frequency = f 
Duty Cycle = D 
RF Choke Inductance = L\ 
Shunt Capacitance= C\ 
Load Resistance = R 
Series Load Capacitance = Co 
Series Load Inductance = Lo 
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S5 










Si 



rkb 3 = ~[m±mR ]dt 
rkb 4 = dt 


Calculation of third set Runge-Kutta constant: 

YYYi = YYi + r -f~ 

YYY 3 = YY 3 + r - fa. 

yyy 4 = yy 4 + r f- 

rkci = |^df 
rkc 3 = 

rkc A = Eg&tfe 

Calculation of fourth set Runge-Kutta constant: 

yyyy 4 = yyy + T -f~ 
yyyy 3 = yyy 3 + 4* 
yyyy 4 = yyy 4 + ^ 


rfcd : 


rkdi = jf-dt 

, 3 = 

rlcd 4 = ?Adt 


Solutions at time t=t+dt by Runge-Kutta fourth order method: 
Y 3 = Y\ + h[rkai + 2rkbi 4- 2 rkc 3 4- rkdi] 


Y 3 


y 3 4- l[rka 3 4- 2 rkb 3 + 2rkc 3 + rkd 3 ] 
y 4 = y 4 + ^[r/ca 4 + 2rA:6 4 4- 2 rkc 4 4- rfcd 4 ] 

= n - n 

t=t+dt 






c ’ 
V J 


PRINT 

t, Yu y 2 , y 3 , y 4 , i 


N=N+1 


(Differential equations for switch-off mode will be applicable) 
Calculation of first set Runge-Kutta constant: 
rkai = 

rka 2 = Ys dt 
rka3 = -[ =Xh±h±XzR ] dt 

r/ca 4 = -Q-dt 


Calculation of second set Runge-Kutta constant: 

YY l = Y x + rJ f^ 

YY 2 = Y 2 + d f z 

yy 3 = y 3 + ^ 
yy 4 = y 4 + ^ 

rkb, = ^Ydt 


rkbi = 
rkb 2 = 


YY^+YYj+YY^R ] 

- ^ dt 

Co al 


Calculation of third set Runge-Kutta constant: 
YYY 1 = YY x + r -^ 

YYY 2 = YY 2 + r -f- 
YYY 3 = YY 3 + r -fL 

yyy 4 = yy 4 + 4* 


vs-yyy> 


v j 
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Ri and ERR2 and ERR$ and ERR4 > LOe 















* PROGRAM FOR 

* USING STATE 


EXACT ANALYSIS 
EQUATIONS . 
double precision 
doub le precision 
double precision 
doub U precision 
open( unit-21, file- 
open(unit=c2, f i 1 9 = 
open ( un i t = £3 , file- 
open ( un i t =84 , f i 1 e = 
open(unit=£5, f i 1 e= 
vs=1 OO . OdO 


OF CLASS-E AMPLIFIER (WITH DIODE) 


y t 1 f q ) , 
rke(^|) ; 
vs,r, 11 
er r ( 4 ) 

= ' 1 1 . i ' 
= ' 1 1 
T t 
1 1 
t T 


y t 2 < 4 ) 
r kb( 4 3 
, 1 O , c 1 


) 


vc ) 
i o ' ) 
vco' 
1 s ' ) 


y(4)>yy(4),yyy(4),yyyy{t$ 

rkc(4),rkd(4) 

C0it ,tt > dt , f ,tpd,.d - 


f=*5G . 0 0 d 3 
d= 0 . 50d0 
11=5. 0d-3 
r s 3 . 08 9 d 0 
c 1 = 0 .2Q3d-6 
1 o = 9 0 . 1 3 d" 6 
co = 0 . I 3d~ 6 
t pd= 1 . OdO/f 
tt-tpd*d 
dt = 1 . Od-9 


1000 
c FOR 


1 0 
100 


m = 0 

K do 1 i = 1,4 
y ( i ) = 0 . 0 d 0 
ytS( i ) = 0 . OdO 
com inue 

t = 0 . 0d0 
n = 3 


S - ON 
do 10 i =1,4 
y t 1 C i ) « y t £ i i ) 
cont i nue 
y ( 2 ) = 0 . OdO 
r k a ( 1 ) = dt * ( vs/1 1 ) 
rka ( 3 ) «-dt » ( y ( 4 1 /lo4*y < 3 >*r/lo ) 
rKa(4)=dt*(y(3)/co) 
yy ( 1 »=y ( 1 )+rka( 11/2, OdO 
yy(3)=y ( 3 ) +r k a ( 3 ) /2 . OdO 
yy (4 )=y (4 )+rka<4>/£. OdO 
r kb ( 1 >«dt*M vs/1 13/* 
rkb(3)=-dt* l y y ( 4 3/3. o + yy(3)*r/lo) 
rkb(4)=dt* C y y ( 3 ) / c o ) 
yyy Cl ) = y y C 1 3+rkbC 13/2. OdO 
yyy(3) = yy ( 3 )+rkb< 3 )/E. 0d0 

yyy{4)*yy(4)+rkb(4)/£,0d0 
, r kc f 1 3 =d t * ( vs/ 1 1 ) 
rkcl3)*-dtMyYyM>/lo+YYY<3J*r/lo) 
rkc( 4)=dt*Cyyy133/co3 
yyyyCI ) — y y y ( 1 )^rkc( 1 ) 
yyyy(3)=yyy(3)'i’rkc(3) 

V y y y C 4 ) =y y y ( 4 1 + r k c C 4 ) 
r kd ( 1 )-dt * ( vs/11 1 

rkd(3}=-dt*(yyyy(4»/lo+yyyy<3>’< ; r/lo) 


.rkd(4)=dt* (yyyy(3)/co) 
y<l)=y(n + (t.OdO/6.0dO)*(rkaU)+£ 
y(3)=y(3)+(t.OdO/6.0dO)*(rka(3)+2 
y(4)=y(4)+( 1 . OdO/6. 0d0)*(rka(4)+2 


. 0d0*rkb( 1 )+£ . 0d0*rkc( 1 )+rkd(1 ) ) 

[od0*rkb(3)+£.0d0*rkc(3)+rkd<3)> 
0d0*rkb(4)+£. 0d0*rkc(4)+rkd(4) ) 


ss=y( 1 )-y( 3) 
t = t + d t 

i f t m . ne . 1 ) go to 2 
if(n.eq.S) then 
n=0 


writ elSl , * H i y C U 

write(22**'3t#yC2l 
wri t e(c3, * ) t , y ( 3) 



°il 


write(24,*)t , y ( 4 ) 
write(25,* : )t,ss 
endi f 

n = n + l » 

2 if(t.le.ct) 90 to 100 
c FOR s _ OFF 

20 0 rks(1) = dT ; Mvs/ll-yC 8 )/li) 

rk a ( 2 ) c dt M y ( t ) / c 1 -y ( 3 ) /c I } 

rk a ( 3 ) =d t * < ( y ( 2 ) -y ( 4 ) )/lo-y(3)*r/lo) 
r k a ( 4 ) = d t * y ( 3 ) / c o 
do 5 1=1,4 

5 yy l l)=yC i J^rkal 1 )/£, OdO 
rkb C 1 )=dtt{vs/ll-yy(E)/in 
rkb(2)=dt*(yy(1)/c1-yyC3)/c1) 

r k b ( 3 ) = d T * C (yy(2)-yy(4) )/lo-yy(3)* : r/lo> 
rkb( 4 )»dt *yy ( 3 )/co 
do 6 i = 1 , 4 

6 yyyU>=yy(i>+rkb<i)/ 2 . 0 d 0 
rkc< 1 )*dx*<vs/] 1 --yyy( 2 >/Il) 
rkc( 2 )*dx*{yyy( t >/ct-~yyy( 3 )/cl > 

r k c ( 3 ) =dt * t i yyy(2)~yyy(4) )/lo-yyy(3)*r/lo) 

rkc(4)=dt*yyy(3>/co 

do 7 i = t , 4 

T yyyy(i)»yyy(i)+rkc(i) 

rkd( 1 )=dt^{\ # s/Il-yyyyCE)/l I ) 

rkd(2)=dt*(yyyy(t)/cl-yyyy(3)/c:1) 

rkd( 3 ) =dt * { { yyyy ( 2 ) - y y y y (4) )/lo~yyyy(3)*r/lo) 

rkd(4)-dtt. yyyy(3)/co 

d 0 8 i « 1 , 4 

3 y ( 1) *y C i ) + ( 1 . OdO/ 6 . OdO ) * t rka< i > + 2 . 0d0*rkb < 1 ) +2 . 0 d 0 *r kc ( i )+rkd( i > ) 

s s = y ( I ) - y ( 3 ) 
i *t + dr 

ir{y(2).it.O OdO 5 then . 

y ( 2 ) « 0 . 0 d 0 
end i f 

1 f { m , n * , 1 ) go to 3 
i f ( n . e q . 2 ) t h e n 
r 1 = 0 

wrixe( 21 ,*)t # y( 1 ) 
wr i t e ( 22 , + > x , y ( 2 ) 
wr 1 1 e ( 23 , » ) t f y t 3 ) v< * 
w p i t e C 2 4 , * ) x , y I 4 ) 
wr i x e ( 25 , * ) x ,ss 

end i f 
n = n+ 1 

3 i f ( t . le.tpd) go to 200 
do 12 is I , 4 

y t 2 { 1 )*y ( i ) 

err(i)*dbs<yx&(i ) - y t t ( i ) > 

12 continue 


i f (err( 1 ) 

QT , 

0 . Id-8) 

go 

to 

1000 

i f ( err ( 2 ) 

gt . 

, 0 . 1 d- 8 ) 

go 

to 

1000 

i f ( err ( 3 ) 

■ g* ■ 

. 0.id-8) 

go 

t O 

1000 

i f ( err t 4 ) 

. gt . 

. 0.ld-3) 

go 

t O 

1000 

in = m+ 1 






i f C m , 1 1 , 2 ) 

go 

to 1000 





c 1 ose ( 2 t ) 
closf (2c ) 
c 1 o s e C 2 3 1 
c 1 05 e ( 24 > 
closet 23 ) 
stop 
end 



Appendix B 

Flow-Chart For Analysis of 
Class-E Amplifier With 
Non-Sinusoidal Output 
Current 


( START ) 



3 

r _ 


CONSTANT 

7T = 3.1415927 

_ 


1 

E . 


READ 

Number of Harmonics in Output Current 
Supply Voltage = Vs 
Operating Frequency = f 
Shunt Capacitance= C\ 

Load Resistance = R 
Series Load Capacitance = Co 
Series Load Inductance = Lo 


NN 


L - 

* 

f 


w = 2irf ; i = 0 

J 

L 
































i=NN+NN+l 
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c 

c 


PROGRAM FOR ANAL V 8 I 
OUTPUT CURRENT . 

parameter ( nn=l () 
doub 1 e preci si on 
double precision 
double preci sio n 
doub 1 e precision 
external fOlaaf 
open (unit=Rl .file 
( uni t 3 c!£ ( f i 1 o 
( uni t =c3 , f i } e 
( un i t =84 , f 3 3 e 


OF CLASS-E AMPLIFIER WITH NON-SI NUSO 1 DAL 

nmax = ( nn+nn+ I j , ia=nmax , ix=nmax , n=nmax ) 
vs, f,c 1 , 1 o , co , pi ,r,t ,tt , t pd , dt 
x ( nn ) , w , vc I , vl , sum! , sumE, q , ql , qq , i 0 
y(n,n),yi(n,nl,b(n,1 I , a ( n , 1 ) ,ctn, 1 ) 
iom(nn) ,phi (nn) 


open 
open 
open 
o p e n 
open 
p i c 3 . 


(uni t=£S, f i 1 e= 
(unit=Z9,fi.le = 
1 4 1 S9E70d0 
vs=1 00 . OdO 
c 1 = 0 . E 0 6 S d- 6 


'ht . 

' y i ' 

' b 1 ' 
'yn 
c 1 ' 


V c 
) 

) 

' ) 
7 


1 


ht . io ' ) 


1 

2 




4 


5 


6 


1 0=367 ; 2 0 E 8 (I d — 6 
co® 0 , 0E84E0 d-6 
r = 2 . 884 d 0 
f=50 . 0 0 d3 


w=£ . OdO* pi *f 


do 1 i = 1 , nn 

x(i>*i*w*io-(l . OdO/ ( i *w*co ) ) 

y t t , 1 )*-n . 0d0/fUi»*8l ) + U . 0 d 0 - ( - 1 ) * * < t ) ) 

do £ i = ( nn+ 1 ) , ( nn+nn ) 

y C i , I } = p i / ( i -nn ) 

y((nn+nn+1 ) , 1 ) = ( p i * *E ) /£ . OdO 

do 3 1=1,00 

do 3 J «e ,, ( not 1 ) 

ifU . eq . < J~1 > ) t hen 

y ( i , j ) = p i * w * c t *»■ r 

else 

y ( i ,))*-(((- 1 ) *M» C i+ J -1 ) >-1 . OdO ) / C i **£- C J-1 ) + *£ ) 

end! f 
continue 

do 4 i » t > nn 

do 4 j = ( nn + 2 ) , < nn+nn+ 1 ) 
if(i . eq . ( j -nn- 1 ) )ihen 

y <-i , j ? * - t p 1 / ( £ . OcfO + i > > +pi *w*c 1 *x < i ) 
clce 

y ( i , j 7 = 0 . OdO 
end i f 
coot i nue 

do 5 i«(nn*J7,(nn + nn7 

do 5 j * 2 , ( nnT 1 ) 

if{< i-nn) .eq. ( j - 1 ) )then 

y ( i , j ) * ( p i / ( £ . OdO * < i-nn) ) ) -p i *w*c1 *x ( i-nn) 


else 

y C t , J ) = Q . 0 d 0 
end 1 f 
cont 1 nue 

do 6 i»(nn+l J, (nn*nn> 
do 6 3*={nn+2), (nn+nn+l ) 
i r t ( 1 -nn > . eq . ( J-n n-1 ) )then 
qq® pifw^cl^r 

y ( i , j ) = qq- < C-J j~nn-l ) ) * < * 

else 

q=* £ < — 1 } *.+ ( J -nn-I ) ) * < 1 . OdO-M - 
ql «=- I . OdO + i r 1 > **l j-nn-l -t-i-nn) 
y(i, 3 )=-qlMi-rm)/U j-nn-l )*C 


OdO-(-l ) * *t ( i — nn ) ) / < (1 - nn ) *r- C j -nn 
)**t i-nn) J 7 /C t j“nn-l 7 7 
( ^ — nn ) * *2 7 ~ ( C J-nn-l )** 2 ) ) J-q 


endif 
cont 1 nue 

do 7 is(nn + nn+Ui<nn + nn+U 

do 7 j = E, ( nnti ) ^ 

y ( i , j ) = U . 0 d 0 - ( - 1 j — 1 
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cio 8 i*( nr^nn+1) , { nn + nn+1) 
do 8 j r ( nn+S ) , ( nn + nn+ 1 ) 

8 y(l,j?=-((-i ) *** ( j -nn-1 ) ) *pi /( j - n n-1 ) 

do 21 l = 1 , ( nn + nn ) 

£1 b ( i , I )=0 . OdO 

b( Inn + mvtl ) , i )«2 , OdOtpit-w + cl tvs 
do 10 1 i " 1 f n 

writei Zc t * ) ( y( i > j ) , j -• 1 , n ) 

101 continue 

w r i t e ( £ 3 , * ) ( b ( i , 1 J , i *» 1 ; n ) 

call f 0 1 a a f ( y , i a . n , y i , t x , a , ifail) 

call mul t { n, n, 1 , y i , b, c) 
write ( 85 , * ) ( c ( i , 1 ) , i «* \ , n ) 
do 1 J I i = 1 ,nn 

iom( i ) *=s qrt C c ( 1 + 1 ,1 )H2 + c(l+nn+1 , \ }+ + 2 ) 

1 f t c( 1 + n n + 1 , t ) .It . 0 . 0 d 0 ) then 

phi ( i J*datan(c(i + 1 , 1 ) /c ( 1 +nn+ 1 , 1 1 ) +pi 

else 

phi (1 )=*datan(c( i + 1 , 1 )/c(i + nn+ 1,1)) 
endi f 

111 continue 

t pd= I . GdO/r 

t.t =• 1 . 0 d 0 / ( 8 . 0 d 0 * f ) 

d t ^ 5 . 0 d ~ 3 

t * 0 . 0 d 0 

£0 in t . ge . 1 1 ) then 

v 1 * c ( 1 i 1 ) * ( w* t -p i ) 

&umj=0 . OdO 
5um£=0 .OdO 
do £33 1=1 i nn 

sum 1 a - t c ( ( i+1 ), 1 )/i ) + s i n ( i*w*t )+sum1 
»uin 2 = ( c ( ( nn + i +1 ) , 1 ) /,i ) * ( cos ( i *w*t ) - ( < - 1) ** i ) ) + sum2 
233 continue 

vd=U . t) d 0 / C w + c 1 ) ) * C v l + 5uml + sumS ) 

u r i t e T2 1 , * ) t , v c 1 

else 

vc 1 = 0 . OdO 
writetcl , * )t , vcl 
end 1 f 

lO-U.Odl} 

a o i :c i , nn 

io^io + i om ( i ) * Us i n ( w*t +phl ( i ) ) 

555 con t i nu c 

WNti'tSDj *)t , io 
t=x+dt 

i r ( t , 1 e . ( tpd+dt ) ? go to 20 

closetcc) 

c 1 os e C 2 3 ) 

c l a s e t 2 4 ) 

c I os e ( 25 ) 

close (cH ) 

stop 

end 

c — * 

subroutine rault iip,aa,bb,cc) 
double precision aa(n,n) ,bb(n, 1 ) ,cc(n, 

, double precision sunt 
do 5 1 = 1 , n 
do 5 j=1 , i ip 

s u m « 0 . 0 
do 6 fc = 1 , m 

6 5um=*sum+ a a t i , k ) *bb t k , j ) 

5 ccC 1 , j 3 *>um 

return 
end 



Appendix C 

Flow-Chart For Accurate 
Design of Class-E Amplifier 


C.l With Vs, Po, Q1 and f Constant 


( START ) 

£ 

CONSTANT 
it = 3.1415927 


READ 

Number of Harmonics in Output Current = N 
Output Power = Po 
Quality Factor = Q1 
Supply Voltage = Vs 
Operating Frequency = f 





104 













































10<J 



PRINT 


Accurate values of Ci, Co, R, Lo 
i ~ ~ — 

( STOP ) 


C.2 With R, Q1 and f Constant 

The flow-chart of this section will be same as the flow-chart of previous 
section, except the last computation block. The calculated values of R and 
Lo, by approximate design method, are retained. Only values of C\ and Co 
are modified. The flow-chart changes, after connection T; it is given below. 







C PROGRAM FOR Accurate Design With Vs, Po, Q1 and f Constant; (Put n > 1 1 " 
parameter(n-ID) 

double precision il , po , vs, f , cl , c t o 1 d , lo* co , pi ,'r 
double precision w, sc, a! n) , bin) , vcnlnl , vcnpsi In) 
double precision 10 ( n ) , z ops i ( n ) , c , 1 oo 1 d , ro 1 d , coo 1 d 
double precision sal ,bbl ,aac 

double p reel u 1 on pb t ( n 1 , 1 om ( n 1 , sum 1 , fuim2 , gum! , gum l 1 , 15 tim 2 £ 

p i "'3 . 14 1 S9,?7 Od 0 

po = 2 0 0 0 . OdO 

ql ttj: 10,0d0 

V5= 1 0 0 . OciO 

f ~5 0 . 0 0d3 

i 1 “pO/V5 

w = 2 . 0 d 0 * p i * f 

phi ( 1 )=datan( -2. OdO/pi )+pi 
c 1 c * p o / ( p i * w * < v s * * 2 ) } 
i oiti ( 1 ) = 1 I /dsinCphi C 1 ) 3 
r=po*2 . 0 d 0/ ( i om ( 1 ) * i om t 1 ) ) 
lo e q] *r/v/ 

ac=l * Od.O/ ( pi *w*ct ) 

a( 1 ) = ac*(2. 0 d 0 * i l+0.5d0*pi*iom( 1 ) ^dcos( phi ( 1 3 ) ) 
sum 1 =2 . Od 0 * i om ( 1 )*dcos(phi ( I ) } 

b( 1 )®&c*(-pi*il-0.5d0*pi*iom( I.)*dsin(phi( 1 3 J-suml 3 
do 1 1=2, n 

c^dfl oat ( i ) 

sum! = ( C c — I . OdO 3 * + ( c-H 3- I . OdO 3 * l om c t J^dsmcphin ) 3 3/Cc**c-1 .OdO 3 
a C i ) =ac* ( i 1 * ( 1 . OdO- ( -1 . OdO ) **c ) / ( ) + suin 1 ) 

1 cont I nue 
do 2 l=£, n 
c=df 1 oat ( 1 3 

sum ) *-( 1 . OdO - C 1 )^c)=Motn( 1 ) *dcosCphH 1 ) ) /c 

sum2=c * ( - 1 . Od 0 + ( -1 . 0 > *•* C c+ 1 . 0} 3 * i om (1 )*dcos(phi < 1 ) 3 / ( c* *2- 1 . 0) 
b ( 1 )»DC*(*-(pi/c)*il+sum2-t-sum1 ) 

2 continue 

20 ( 1J*sqrt ( a( 1 3**2+b( 1 3 **2 3 / C 1 om Cl 3 3 
co^l . Od 0/ ( w*w* 1 o-w* ( sqrt ( zo ( 1 3 *zo ( 1 3 -r*r ) 3 3 
.writet^,*} 'Approximate values of c I , co , r , 1 o= ' > c 1 , co , r , 1 o 
11 00 do 3 i*1 ,n 
c=df 1 oat C i ) 

vend ) = sqrt U( 1 ) **E + b C 1 3**2) 

20 ( 1 Jasqrt (rMt2+(c=i:wdo-( 1 . OdO/(c*w*co> > >**£) 
zopsl ( 1 3=atanC <c*w*lo-t 1 . OdO/'lc + w^co) ) )/r) 
i om C 1 3 =vcn C i ) /to <’i ) 

1 f C b t i ) . 1 1 . 0 . OdO ) then 

venps i t i .3 *at an ( a ( i )/b( i ) 3 • p i 

else 

vcnpsH 1 3 = at an ( a ( 1 3 / b ( 1 3 3 
end! f 

phi ( i )*vcnpsl ( i ) -zops 1 ( 1 3 

3 cont 1 nue 
sum 1=0 .OdO 
sum 2=0 . OdO 
do A 1 =2 , n 

c-df 1 o at (13 . _ t - 

sum 1 =sum l + ( l . OdO-< - 1 . OdO ) ** i ) * 1 om< i J *dcos ( pm 1 1 ) ) /c 

sumH = sum2+ i om ( i )*<jsin(phi ( i J ^ 

A cont inue 

aal = — 0 ,5d 0* ( pi*11 +sum 1 3 
b b 1 = i 1 - s u m 2 

lomt 1 ) = sqrt (aal + * : *Z1 

if ( aa I . It .O.OdO ) t hen 
phi Cl ) -at antbbl /aal 3+pi 
else 

■ phi ( 1 3 =at an( bbl /aal 1 
endif 
cl o 1 d=c 1 



6 


S 


9 


13 


coo 1 d= co 

rold=r 

loold=*lo 

I ,0d0/(2.Qd0*pi*w*vsJ 
sum ! - 0 . OdO 
£Ulfic = 0 . 0 d 0 
do G i « 1 , n 
c~d NoatliJ 

5uml*suinl + u . 0 cJ 0 — C — 1 . 0 d 0 7 * * i ) * i o m ( t ) * d 5 I n { p h i C i 3 
s urn c = s u m2 + ( (-1 . 0 d 0 7 * * i 3 * p l * 1 o m ( i )^-dcos(phi C i ) 7 / c 
cent i nu? 

c 1 -• C 0 . 3 d 0 t. p i * p i 4 : 1 1 + s u tn 1 - s u m £ 7 * a a c 
ac**l . OdO / ( pi nw+-cl ) 
do 9 i 21 1 > n 
c=df lost < i ) 


s u m 1 30 0 . OdO 
sum 11 =0 . OdO 
sumEE™ 0 .OdO 
d o 3 J = 1 , n 
e*=cif 1 o at ( j ) 
ifti.eq.j) then 

suml =5um) +0 .5d0*pi*iom( j )*dcos(phi ( j ) ) / c 
sum! I =sum 1 1-D.5d0*pi*iom( j )*dsinlphi ( j ) 5/c 

sumEE^sum ££+((-1.07* *17*11 . 0- ( - 1 . 0)**i)*iom(j)*dco5(phi(j))/<c*e) 
else’ 


sumE^c^E-e* *E 

. OclGi^^C )~i 

5um3» , £*Mc + *E-e* + E) 


OdO i’i'ioiM 


3 ;*dsimpni \ J 


/sufnc 


sumi l=sunrn+c*((-i .o)=^Mi + ])-l . 07*lom(j}*dco5(phi(j7}/sum3 
5 um££= 3 5 Uin£E+ ( (-1 . 0 ) *•+* 7 * C 1 .0-1-1 .07**1 ) * i om ( j > ^dcos ( ph i ( J ))/(c*e) 


end i f 
cont inue 

a ( i ) = a c * ( i I * ( 1 . 0 d 0 — ( ~ 1 . 0 d 0 7 * * i )/(c**2)+sum! 7 
b ( 1 J = ac ! M ( — p 1 / c 3 * 1 1 + Bum 1 I +5 uid2£ J 


continue 

20 ( 7 7 =sqrt ( * ( 1 7 * E + l> ( i 7**c)/iom( 1 7 

sum 1 - 0 . OdO 

do 73 1=1 , n 

sum 1 =sum 1 +.1 om ( i 7 **£ 

cont lnue 

r=£ . OdO*po/sum 1 

1 o**q 1 *r/w 

.co^l ;■ OdO/ ( w* w* 1 o-w* J C sqrt ( zo C 1 )*zo(l ) ~r*r > ) > 
if C C'abs (cl-cl old ))• ot . 0 . 000001 d-6 ) go to 1100 
l fCCabcC co-coold )). gt . o . ooooo l d-6 ) go to iioo 
i f ( Cabs ( 1 o- loold ) ) . gt . 0 . 0 0 0 0 0 1 d-5 7 go to 1100 
1 f ( ( abs ( r-ro 1 d ) ) . gt . 0 . 0 0 0 0 0 1 d- 6 ) go to 1100 
wr i T e ( * / * ) "Accurate values of cl ; co , r , 1 o =',c1,co,r,lo 
stop ' 

end 



c PROGRAM FOR Accurate Design With R, Qi and f Constant; 


(Put n > 1 ) 


paramet &r(n=!0) 
double precision 
double precision 
double precision 
double precision 
double precision 
p i - 3 . 1 A 1 3 9 2 T 0 d 0 
pos? 000 . 0 d 0 
q 1 * 4 0 . 0 d 0 
v s = I 00. Odd 


i \ »pO; vs , f i q 1 , c 1 , lo, co,pi / r 
w,ac,a(n) ,b(n),vcn(n) , vcnpsi (n) 
xo(n) , zopsi (n) , c , coold f clold 
a b 1 / b h 1 , a a c 

phi (n) , i o m ( n ) , s u m 1 , ? u r 8 , sum 3, sum] 1 , sumH? 


f = 50 . 0 0 d 3 


i 1 s po/vs 
w^£ . 0d0*p i * f 

phi i 1 ) =dat an ( ~2 . OdO/pi )+pi 
c 1 =p o / ( p i * w * ( v s * * 2 ) ) 

.iom( 1 } = i l/dsin(phi( 1 ) ) 

r = po*c. 0 d 0 / ( i o m ( 1 )*iom( 1 ) ) 

lo*ql*r/w 

ac= 1 . OdO/(pi*w*dJ 

a M )=ac* ( £ . OdO* i 1 +0 . 5d0*pi *iom( 1 ) *dcos< ph i ( 1 3) ) 
suni1~2.0d0*iom( 1 ) *dcos ( phi ( 1 ) ) 

b( 1 ) =a c * ( ~p i * i 1 - 0 . 5 d 0 •+ p i * i oro ( 1 }*dsin(phi( 1 ) ) - s u m 1 ) 
do 1 i =2 , n 
c = df 1 oat ( i ) 

sum i *= c { c - 1 . ocio ) ** ( c+ 1 ) - 1 . ocio ) * lom c i Hdsmtphit i > ) . odo) 

a ( : i ) =ac* ( i 1 * C 1 .0d0-C-l . 0d0)**c)/Cc**2) + sutnl ) 

1 cont Inue 

do 2 1=2, n 
c=df 1 oat ( i ) 

■5urr»]=-( ] . 0d0~(-l )**c)*iom( 1 )*dcos(phi ( 1 ) )/c 

5ur»?=c : ^ C - 1 . 0d 0+ ( — 1 . 0 ) ** ( c + 1 . 0 ) ) * i om ( \ ) *dcos ( ph i ( 1 ) ) / ( c* *2-1 . 0 ) 
b(.i J»ac*<-tpi/c)*i 1 +sum2 + suml 1 
2 cont inue 

zo n >=sqrt < a! 1 )**E+b(l 1 **2 ) / < i om ( 1 M 

co-1 . 0d0/(w*w+lo~w*(sqrt<zo( l ) *zo ( I ) ~r*r > ? ) 

write! * , * ) ' Approximate Values of cl ,co,r,lo--' ,d ,co,r, lo 
1100 do 3 1 = 1 , n 
c=df lqat ( i ) 

vend )=sqrt (a( i )*^ : c+b( 1 )*< : c) 

loi i )=sqr t(r**2*(c*w*lo~( 1 . Od 0/ ( c*w*co ) ) ) *• 2 ) 

z ops i ( i > =at an ( < c*w* 1 o- ( 1 . 0d0 / ( c^wt-co ) ) ) /r ) 
lomd )*vcnC i )/ 2 o(l ) 
trcbcii .it .o.odoithen 
venpai ! 1 )*5tan( a( 1 )/bd } )+pi 


else 

venps 1 ( i )=atanta( i ) /b ( 1 ) ) 
endif 

ph i C i ) » venps i C i)-topsl(i) 

3 continue 

surnl = 0 . 0d0 

sum2=0 . Od 0 
do A 1 =2 , n 
c=d f 1 oat < i ) 

sum l^sutn 1 +< 1 . 0 d 0- ( -1 . Od 0 ) ** i ) * 1 om ( i 1 +dcos C phi < i >> /c 
5 um 2 =sum£+ i om < i ) *ds 1 n < ph i ( i ) ) 

A cont inue 

aa 1 *-0 * 5ti0* ( p 1 * i 1 tsurol ) 
bb 1 = i 1 -sutnS 

ioiri( 1 ) -sqrt ( a a 1 **2+bb 1 **H > 
i f ( aal . It . 0 . OdO } then 
phi t U = atan!bbl/aal )+pi 
else 

phi ( 1 1 =at an ( bb 1 /aal 1 
end i f 
cl old-cl 
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coo 1 d = c o 

aac=l . 0d0/(cf.0d0*pi*w*vs) 

suml» 0 . 0 d 0 

s U rr< £ = 0 . OdO 
cl o 6 1 = 1 , n 
c*df 1 o a t ( i ) 

ii u iii 1 * s u in 1 -i l 1 . 0 d 0 - ( - 1 . o d 0 M s * 1 ) * l o m ( 1 ) * d s l n ( p h i c i } ? / ( c * * £ 3 
suiric = suii\2+( (-1 . OdO ) * *• 1 ) rp 1 * i om C i )*dcos(phi C i ) )/c 
6 continue 

ci-u?.5dOTpi%pi*n +sufn l - sum c?*fi<ac 

ac=l . 0d0/(pi**w*c1 ) 

do 9 1 = 1 . n 

c=df ] oatllJ 

s u rrt 1 -• 0 . 0 d Cl 

sum M =0 . OdO 

sumcic = 0 .OdO 

do S j = t , n 

e = dr 1 oat ( J 1 

i f ( 1 . e q . 1 ) then 

sum 1 = suml +0 . 5d0*pi *i om( j ) *dco 5 ( phi C j ) )/c 
sum 11 =sum II -0 . 5d0*pi *iom( j ) *dsin( phi C j ) )/c 

sum22 = sun»22+ ( ( -1 . 0 ) **i )*( 1 . 0-( -1 . 0 ) ** i ) * i orn ( j ) *dcos ( ph i ( j ) ) / ( c*e ) 
else 

surn£ = c**2-e**£ 

5uml*5uniM-( ( **■ 1 . 0 d 0 ) * * ( 1 + j ) - 1 . OdO ) -+iom( j )*dsin(phi(j) )/suri 2 
sum3=e*( c+*£-e**£> 

sumn = sum 1 1 +c* ( ( -1 . 0 J + j 1-1 . 0 ) *1 om t j J *dcos ( phi ( j ) ) / sum 3 
« umS£-2.uni£2+ ( (-1 . 0 ) * * e ) * ( 1 . 0- ( ~ 1 . 0 ) ** i ) + lom( j )+dco3(phi ( j ) ) / ( c*e ) 
end i f 

8 continue 

M i )=ac*( i i *( 1 . OdO-M-J . 0 d 0 ) ** i ) / ( c**E ) +sum I ) 
b ( i ) = ac* ( ( ~p 1 /c J * 1 I tsuml 1 + sum2c! ) 

9 continue 

20 ( 1 ) =sqrt ( a ( 1 }**2 + b( I ) * *2 ) / 1 om ( 1 ) 

c o = 1 . 0d0/(w*w*lo-w*(sqrt (io( 1 ) *zo ( 1 } - r * r ) ) ) 

sum 1=0. OdO 

do 1 3 i= 1 , n 

sum 1 =sum 1 + 1 orn ( i ) **£ 

13 continue . 

po s n ; 5Utnl /5 . OdO 
vs^po/i 1 

i f (( abs ( c 1 -cl old )) .g t . 0 . 0000001 d-6 ) go to 1100 
i f ( ( abs ( co — coo 1 d ) ) \ g t . 0 . 0 0 0 0 0 0 1 d — 6 ) go to 1100 

urlteC*,^) 'Accurate values of c 1 , co , r , 1 o , vs = ' , c 1 , co , r , 1 o 
stop 
end 
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Flow-Chart For Exact Analysis 
of Class-E Push-Pull Amplifier 
With Two Modes of Operation 
(Using State-Equations) 

( START ) 

X 

READ 

Supply Voltage = Vs 
Operating Frequency = 

Duty Cycle = D 

RF Choke Inductance = Li = L 2 
Shunt Capacitance = Ci = C 2 
Load Resistance = R 
' Series Load Capacitance = Co 
Series Load Inductance = Lo 

J_ 

tpd = Time-period = j 
tl = Time duration for mode-1 = j 
t2 = Time duration for mode-2 = 
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program for push-pull class -l amplifier cuith d*o.5). 

double precision y(6),yy(63 / yyy(6),yyyy({S) 
double precision rkc(&) , r k b ( 6 > , r k c ( 6 3 ; rkd(6 J 
double precision vs, r, 11 , 12^1 lo ; co/c , f ,tpd 

double precision t 1 , r 2 , s si , s s 8 , d 
double precision v t 1 C 6 3 , y t 2 < 6 ) , err ( 6 ) 
p e n i u n i i = 2 1 , f i 1 e - 7 1 t . i 1 ; 3 
open ( uni t -22, f i 1 s = ' i r . vc t r ) 
open(unit=23;rile='tx . i 2 ' ) 
open(uni t - 2 4 , f i i e - ' t c .vcc" ) 
o p e* n ( u n i f i 1 c- = ' c t . io' ) 

open(unit=t:*6, f i 1 e = ' t t . vco ' ) 
ope n (unit-27, file-' tr . i s 1 ' ) 
open ( un i t -28 , f i 1 e = ' x i . i sc* 1 ) 
v s=1 00 . OdO 
r-S . 7 6 8 d 0 
11*5: 0d-3 
18-5 . 0d-3 
cl -0 . 20264d-6 
c2=0 . 20264d-6 
lo-183 .60! d-6 
co = 0 . 0 62 37 4d- 6 
f = SO. 0 d 3 
t pd= 1 . OdO/f 
d-O.SOdO 
1 1 * t p d * d 

x2~t pd* ( 1 . OdO -d ) 
do t i » 1 , 6 
y ( i ) = 0 . OdO 
.yt2(i)*0.0d0 
I continue 

d t = 1 . Od-9 

JD s 0 

nn = 5 0 

1 0 0 0 t = 0 . OdO 

do 10 i s l/6 
yx I ( i‘) = y t£ ( 1 ) 

1 0 cont i nue 

FOR CASE OF SI _ OFF * S2 - ON 


100 y ( 4 ) - 0 . OdO 

r k a C 1 ) s di:^( vs/ ( J1 -y(c)/l 1 3 
, rka(2)=dt , * e (y(1}’/cl-y(5)/cl)- 
rka(3)=dt : Mvs/12) 

rkaCS)*dt*( (y(23"**y(63 ) / 1 o — y(5)*r/lo) 

rka(6)=dT*y(5)/co 

do 5 i - 1 , 6 

5 yy( i ) = y ( i ) + rka( i ) / 2 . OdO 
rkb ( 1 ) c dt * ( vs/ 1 t ~y y ( £) /l 1 ) 
rkb(2)=df*(yy C 1 3/cJ-yy (5)/cl ) 
rkb(3)«dt*< vs/12) 

rkb ( 5 ) =dt * ( {yy(2)-yy(6))/lo-yy(5)*r/lo} 
rkb { 6 ) =d t * y y ( 5 ) / co 
do 6 i = 1 i 6 

6 yyy( i )=yy ( i )+rkb( i )/2 . OdO 

rk c ( 1 3 «dt * ( vs/ 1 1 ~y y y l 2 ) / 1 1 ) 
rkc ( 2 )=dt * ( yyy < 1 ) /c 1 -y y y ( 5 3 /c 1 ) 
rkcOi-dt + Cvs/lE) 

rkc ( 5 3 = dt * ( ( yy y ( 2 ) ~y y y ( 6 3 3 / 1 o-y y y ( 5) *r/l o 3 
rkc C 6 )*dt»yyy(5)/co 
do 7 i- l , 6 . 

7 yyyy'(i)®yyy(i)*rkc( i ) 

rkd( t )*dt* ( vs/ It -yyy y (£>/l 1 ) 
rkd( 2 ) *di * < yyy y U ) /cl -yy y y C5 ) /cl 3 
rk d ( 3 ) =di * ( vs/1 c ) 



r k d { 5 ) - d L * ( { y y y y { £ ) - y y y y ( 6 ) )/}o-yyyy(5)*r/lo) 

rkd(6}=dt*yyyy(5)/co 

do 8 i = 1 , 6 

6 y < i >=V < i > + M . OdO/6 OdO)Mrka( i M2 . OdO*rkb< i >+£.'. OdO*rkc< i >+rkdC l 
I = r + d r 

i f i y C 2 ) . 1 & . 0 . OdO f ti'irii 
• y ( 2 ) -= 0 . 0 d 0 
end i f 

*si~y( 1 ) - y ( 3 ) 

■ ss£*y (3)+-y C£) 
y ( 4 ) = 0 . 0 d 0 
i f C m . n e . I ) go to 3 
i f ( nn . eq . 50 i then 
writ e < 21 , * ) t , y ( 1 > 
w r i t e ( 2 2 , f ) t , y ( £ ) 
writeC23,*)t,y(3) 
ur i t r ( 24 , * ) t , y ( A ) 
uri t * ( 25 , * ) t , y ( 5 ) 
w r i t e ( 2 6 , * ) t , y ( 6 ) 
wr i t e ( 2 7 , * ) t , ss 1 
wr i r & ( 28 , * ) t , ss2 
nn- 0 
. end i f 
nn=nnf 1 

3 if (*c . It . t 1 ) go t o t 0 0 

c FOR CASE OF SI ON $ Sc* - OFF 


300 rka{})~dt*(vs/li) 

rka(3)«dt*<-y(4)/I 2+vs / 1 2 ) 
rka(4)»dt*< y(3)+y (5) )/c2 
r k a ( 5 ) = - d t ( y ( 4 ) + y ( 6 > ) / 1 o •+* y ( 5 ) : * r / 1 c ) 
r ka ( 6 ) «dt *y ( S ) /co 
do 9 i=1 , 6 

9 y y ( i ) = y ( i ) + r k -i ( i ) / £ . 0 d 0 

r k b ( 1 )«dl*( vs/1 I ) 
rkbC3)=di*(-yy(4)/12+vs/I2) 
rkb ( 4 ) =dt + ( y y ( 3 ) + y y ( 5 ) ) /c2 

rkb(5) 3 =~dt*((yyM) + yy{6)>/lo-*-yy(S)*r/lo) 

rkb(6)=dt*yy(5)/co 

do 6 6 i - \ j 8 

es . yyyCi)*yy(i)+rkbU)/2‘. OdO 
rkc( 1 )-dt* ( vs/1 J > 
rk c ( 3 ) *d t * ( -yy y {4 . ) / 1 £+vs/12 > 
r k C( 4 )’=dt » ( y y y ( 3 ) +-y y y ( 5 ) ) /c2 

r k c ( 5 ) = r d t * ( { yyy(4)+yyy(6 ) )/lo+yyy(5)*r/lo) 
rkcCSl-dt^yyyl 51/co 
do 77 i* 1 , 6 

77 yyyy(i)=yyy(i)+rkc(i> 

rk d ( 1 )«dt*(vs/l 1 ) 
rkd(3)=dt*(-yyyy(4)/12+*vs/12) 
r k d ( 4 ) =•* d t * ( yyyy(3) + yyyy(5) ) / c 2 

rkd(5)--dt* : ( (yyyy(4) + yyyy(6) )/lo + yyyyC5)*r/lo) 
rkd{&)= : dt*yyyy(5)/co 

do 88 i=1 ,6 ^ 

S? y ( i )=y ( i ) -M 1 . OdO/ 6 . OdO )* ( r k a( i H 2 . 0-d0*rkb ( i ) +c . 0d0*rkc ( i ) +rk d ( i] 

y(2)=0.0d0 
t +dt 

if ( y (4 ) . 1* . 0 . OdO > s th$n 
y ( 4 ) ~ 0 . OdO 
*ndi f 

ss.i *=y ( 1 ) -y ( S ) 
ss2 = y ( 3 ) +y ( S > 
i f ( m - n e . 1 ) go to 33 
i f l nn . eq . 50 ) then 
w r i x # ( 2 1 ,*♦■■) t * y ( 1 1 
wriTe(2£/*)T,y(2J 



33 
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writc-CSl, 

•Hi , 

y (3) 





ur i t e ( £4 , 

* ) t , 

y { 4 ) 





ur i t e ( 2S , 

* > t , 

y (5) 





w r i i * ( c 6 , 

* ) T , 

y ( CO 





ur i t e ( ;? ? . 

; mt, 

ss 1 





wrj v e ( 2 8 , 

* u , 

ss2 





r.n^O 







~ n d i f 







nn-nn-t ; 







i f ( T . 1 6 . ( 

t r+t £ i ) 

go to 

300 



do 12 i = t 

,6 






y t 2 1 i ) = y ( 

i ) 






; j rr( i ) - a b s ( y t 

2 ( i ) 

-ytl ( 

i ) ) 



continue 







i f ( e rr ( 1 ) 

• gt 

. 0 . 

1 d-3) 

go 

t 0 

1000 

L f l e r r ( £ ) 

- gt 

. 0 . 

1 d-S> 

go 

t 0 

1000 

i f l err ( 3 ) 

- gt 

. 0 . 

1 d-8 ) 

go 

t 0 

1000 

l f 1 e r r ( A ) 

• gt 

. 0 . 

1 d-8 ) 

go 

t 0 

1 000 

i fUrrt 5) 

* gt 

. 0. 

1 d-8) 

go 

t 0 

1000 

i f ( e r r ( 6 ) 

• gt 

. 0 . 

1 d-8) 

go 

t o 

1000 


ttt - ill + 1 

i f ( in . 1 t . «? > qo to I 0 0 0 

close(8l ) 

c.los*( £<? ) 

c l o s e* ( c 3 ) 

c 1 o s e ( 2 4 ) 

cl o s e ( •: S j 

close (56 ) 

close ( 2? ) 

c 1 o s e ( c h ) 

st op 

* n d 
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Flow-Chart For Exact Analysis 
of Class-E Push-Pull Amplifier 
With Three Modes of 
Operation (Using 
State-Equations) 

( START j 

z 

READ 

Supply Voltage = Vs 
Operating Frequency = f 
Duty Cycle = D 

RF Choke Inductance = L\ = L 2 
Shunt Capacitance = C\ = C2 
Load Resistance = R 
Series Load Capacitance = Co 
Series Load Inductance = Lo 

! k 

tpd = Time-period = j 





1 ii 




V J 


tl = Time duration for mode-1 = Dd 
= t4 = Time duration for mode-2 = 
t3 = tl = Time duration for mode-3 
Time step for calculation = dt = 

400C 

m=0;N=10;i=0 


i=i+l 


Yi = 0 . 0 ; YT2i = 0.0 


Wo / 18 

— — <T i = 6 


t=0.0;i=0 


i — i+l 


YTh = YT2i 


i = 6 


Y 4 = 0.0 

Solve Following Differential Equations For Mode-1: Si - OFF & S2 - ON 

(Here, Y x = h]Y 2 = V a ; Y z = J 2 ; Y 4 = V&; Y s = I 0 and Y 6 = Vco .) 

= Ys. — Yl . - 

dt L\ L\ 


r > 

3 

v y 





























U7 





PRINT 

h *1, f 2 , y 3: r 4; y 5 

y& 7 1 si, hi 



y 2 = o.o y 4 = o.o 

Solve Following Differential Equations For Mode-2: Si - ON k S2 - ON 
(Here, Y 1 = I 1 -,Y 2 = V cl ;Y 3 = I 2 ; 7 4 = V; 2 ; Y 5 = I 0 and Y & = Vco .) 

dVi _ 7, 
dt Li 

m. - Ji 

dt L/ 2 

dt Lq L 0 

dYe _ Yk 
dt C 0 

(The above equations are solved at time t=t+dt by 
Runge-Kutta fourth order method) 

i s i = y x - y s; j s2 = y 3 + y 5 ; t=t+dt 
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C PROG'S AH FOR CLASS-E TUNED AMPLIFIES WITH 3-MODE OPERATION (D = 2/3). 

double precision y ( G ) ,'y y ( 6 ) , y y y ( 6 ) , y y y y { 6 ) 

double precision rke.( & ) , rkb ( 6 ) , rkc ( 6 ) , rkd ( 6 ) 

double precision vs, r , 1 1 , 12, ci , c£, io, co, t , f , tpd, d 

double precision t 1 , i £, t 3 , t 4 , ss 1 , ss£ 

double precision yt1(&),yt£<6),err(6) 

open ( uni t = £ 1 f i 1 e - ' vn . i I ' ) 

open ( un i t -££ , fi 1 e = 't n . vc 1 * ) 

0 p e n ( u n i l = £ 3 , f i. 1 e = ' 1 n . 12' } 
open(unit-£4,file='tn.vc£ r ) 
openCunit-ES, fi 1 e ~ ' tn. i o ' ) 

Mpen(uniT=£6, f il e = ' rn.vco' ) 
open(unit=£7, f ile*‘tn. i si ' ) 
openCunit*£8, f ile='tn. i s 2 ' ) 
vs- 1 00 . OdO 

f “SO . 0 d 3 
d = 2 . 0d0/3 . OdO 
11*5. 0d~3 
■ 12=11 

cl =0 . 031 6Sd-6 
c2 = 0 . 031 6Sd-6 
r = 1 1.971 OdO ' 
co=0 . 0£839d-6 ' 

1 o=3-S I . 0 65£d-6 
t p d = 1 . 0 d 0 / f 

t 1 = t p d *• ( 1 . 0 d 0 - d ) 
t3 = t’l 

t 2 - t p d / 6 OdO 
t 2 = t p d * ( 1 . 0 d 0 d ) ** 0 . 5 d 0 
t h = t £ 

d t * t D d / 4 0 0 0 . 0 if 0 
iti ~ i ) 
n = 1 0 

do 1 i = 1 , 6 
y ( i ) = 0 . 0 d 0 
yt£ ( i )=0 . OdO 
1 continue 

1 00 0 . t = 0 . OdO 

do 1 0 i«1 , G 

y M ( i ) •= y t: 2 C i > 

10 continue 

i: FOR CASE OF SI - 0 rf" S2 ~ ON 


1 00 y ( 4 )~0 . OdO 

rkaHl^dTt-Cvs/ll-ylSI/ll) 
rka(2)-dt*(y(l )/c1-y(5)/c1 ) 
r k a ( 3 ) ~ d t * ( v s / 1 2 ) 

rka(S)=dt4( ( y ( £ ) - y ( 6 ) )/lo-y(5)*r/lo) 

. rkB(6)-=dt + y(5)/co 
do 5 i 11 1 > 6 

5 yyt i > — y ( i ) + r k a < l ) / £ . 0 d 0 
rkb( 1 )~dt+( vs/li-yy ) 

rkb ( 2 ) = dt * < y y ( 1 ) /c 1 - yy ( 5.) / c 1 ) 
fk b ( 3 ) =dt * ( vs/ 1 2 ) 

r k b ( 5 ) =dt * < ( y y ( 2 ) -y y f 6 ) ) / 1 o- y y ( 5 ) *r/ 1 o ) 
rkb(6)=dt*yy (5)/co 
ci o £> i = 1 , 6 

6 yyy(i)=cyy(iJ+ r kbii)/2,0d0 

rkc( I )»dt*(v€/]l-yyy<2>/ll ) 

r k c, ( £ ) = d t * ( y y y ( 1 ) / c 1 - y y y ( 5 > / c I ) 
rkc(3)=dt*( vs/12) * 

r k c ( 5 ) = d r * ( (yyyl£)~yyy^6J )/lo*-yyy(5)*r/lo) 

rkc(6) x dt+yyy £&)/co 

ao 7 l - 1 / 6 



ilo 


( yy-yy ( i >*yyy ( i ) +rkc c i ) 

rkd( 1 )=dtM vs/] t^yyyy(8)/l 1 ) 
rkd(E)»dt*Cyyyy< t )/c1-yyyy<S)/d ) 
rkd ( 3 ) = dt * ( vs/ i£ ) 

r k d ( 5 ) * d r * ( C yyyyCEI-yyyyCG) )/Io-yyyy(5)*r/lo) 
r k d ( G ) = d t + y y y y { 5 3 / c o 
do 8 i » 1 , 6 

F ' Y ( i ) *= y ( 1 ) + ( l .DdO/6.QdG)Mrka(i)+£..OdO*rfcb(i)+E OdD*rl;c( i ) + rf c!( i } 

} =t +dt 

i f ( y ( 2 ) . 1 t . 0 . OdO ) then 
Y ( d ) - 0 . 0 d 0 
*? n d i f 

«l=y< l )-y(5) 

3s£=y (3)+y (S) 
y ( 4 ) =0 , OdO 
i f ( m*. n ? . 1 ) go to 3 
i f ( n . e q . 1 0 ) then 
n- 0 

wr i t e ( E 1 , * ) t , y ( 1 ) 
writ e (££, *) t , yt£.) 
writ e(£3, *)t , y (3) 
w r i t e ( £ 4 , * ) t , y ( 4 ) 
wr x t * < £S , * > t , y ( 5 ) 
wr i t e ( £ f) j -t : ) t ; y ( 6 ) 
wr ) tei c7; * ) t , ss 1 
wr i t e ( 2 8 * ) t , S s£ 
end! f 
n=n-H 

3 i f ( t .. 3 t . t I ) go to IDO 

-FOR CASE OF SI ON & S3 - ON 


y ( c‘ ) « 0 . 0 d 0 
V ( 4 ) = 0 . 0 d 0 

300 rka ( 1 ) -dt * ( vs/1 1 ) 
rk a ( 3 ) - dt * ( vs / I'd ) 
r k a ( 5 J. - - d t * C y ( 6 ) / 3 o + y ( 5 ) * r / 1 o ) 
r k a ( G ) ~ d tt- { y l S ) / co ) 
yy U ) e y ( 1 )+rk^( I )/£ OdO 
y y ( 3 J = y ( 3 ) + r k a ( 3 ) / 8 . 0 d 0 
yy (5.) = y(5>+rka(5>/£. OdO 

V y(6)=y(6)+rkaC6)/3.0dO 

rkb{1)=dtMvs/lU 

rk h ( 3 ) =dt * ( vs/ v l 3 ) 

r k b ( 5 ) d t * ( y y l 6 ) / 1 o+ y y C 5 1 * r/ 1 o ) 

r k b ( 6 ) = d t * { y y ( S ) / c o ) 

yyy( l ) -y y ( 1 ) +r kb ( I ) /Z . OdO 

y y y ( 3 ) = y y ( 3 )+rkb(3)/E.0d0 

yyy(5)=yy(5) + rkb(5)/‘£.0d0 

y y y { 6 ) - y y ( 6 )**-rkb(6)/£‘-0d0 

rkc(l) x dtMvs/ll> 

rk c ( 3 ) ~d 1 + ( vs/ 1 c ) 

rkc(5)=Tdt*(yyy(6)/lo+yyy(5)*r/lo> 
r k c ( 6 ) « d 1 1- ( y y y ( 5 ) / c o ) 
y y y y ( 1 ) * y y y ( I > **• r k c M ) 
yyyyC3) 5 * r yyy(3)^rkcf3) 
y y y y C 5 ) * y y y ( 5 ) + r k c 1 5 ) 
y y y y ( 6 ) = y y y ( 6 ) + r k c ( 6 ) 
rkd ( 1) -dt + ( vs/1 1) 
rkd(3)=dt*(v£/13)- 

rkdC 5 ) *-dt *■ { yy yy f b )/ 3 yyyy (S ) *r/3 o ) 

r k d < 6 ) « d t * ( y y y y ( S ) / co) 

" y ( j >. sy (.| ) + ( \ . OdO/6 . OdO )*( rka( 1 ) -*2 . OdOtrkb( 1 )-*E.OdO*rkc< 1 )*+rkd(1 ) 

y ( 3 ) =y ( 3 ) + { 1 . OdO/6 . OdO ) * ( rka ( 3 ) +3 . OdO ^rkb ( 3 )+£ . OdO + rkc ( 3 ) •* rkd < 3 ) 

v C 5 ) =y ( S ■) + ( 1 . OdO/6. 0d0)*(rka<5) + 2 0 d 0 + r k b ( 5 ) ^3 . 0 dO < r k c ( 5 ) +r k d ( 5 ) 

y(6) »y(6)*+( I OdO/G. 0cJ0)^(rk*(G) >3, OdO-* rkb( 6 )•*■£. OdO 4 r kc( 6 )•« rkd( 6 ) 

1 «t **-d t 



'SS 1 ay ( { ) *~ y ( S ) 

3s£ = y(SMy(3} 
if (m.ne . 1 ) go to Z 
i f ( n . eq . 1 0 ) then 
n= U 

write(21,*)t,y{l) 

-r i r e l 22 , + ) i; , y ( £ ) 
wpi t e( £3j + ) t , y (3 ) 
wnte(2h*)t,y(4) 
wri t e( £5 , * ) t , y C S ) 
write(26,*)t,y ( 6 ) 
wr i t e ( 27 , * ) t- , s s 1 
wr i t e ( 28 r * ) t t ss2 
end t f 
n « n + 1 

£ i f ( t . ] t . ( t 1 + t 2 ) ) go to. £00 

c FOR CASE OF SI _ ON & S£ - OFF 


300 y < 2 ) « 0 . 0 d 0 

rka( T ) ~ d t * ( v s / 1 t ) 

rk a ( 3 ) -d t * ( ~y M ) / 1 2 + vs/ 1 £ ) 

r k a ( 4 ) = d t * ( y ( 3 ) + y ( 5 ) > / c2 

rka(S)=-dT*((y(4)+y(6))/lo+y(5)*r/lo) 

rka(b)-dt*y(S)/co 

do 9 i *= 1 , 6 

* yy(i)*y(i)*+Tka(i)/2.Qd0 

rkb ( 1 ) = dt * ( vs/1 1 ) 

• r k b ( 3 ) **d t * ( - y y ( 4 ) / 1 2+ v s / 1 £ ) 

»*• k b ( 4 ) ~ d t * ( y y ( 3 ) + y y C 5 ) ) / c 2 
r*kb(5)“-c!t*( (yy(4)+yy<6) )/lo+yy(5)*r/lo) 
rkb ( b ) =dt *y y ( 5 ) /co 
do 6 6 i » ! , 6 

66 y y y ( i ) =s y y ( i ) 4 r k b ( i ) / £ . 0 d 0 

r k c ( 1 ) « d t * ( v s / 1 1 ) \ 

r k c ( 3 ) = d t * ( - y y y ( 4 ) / I Z + v s / 1 2 ) 
r k c ( 4 ) = d t * ( y y y ( 3 ) + y y y { 5 > ) / cZ 

rkc(5)=-dT : *((yyyM)+yyy(6)}/lo + yyy(5)*r/lo) 

rkctGJssdt^yyyCSl/co 

do 77 i - 1 t 6 

7*7 y yyy ( i ) *yy y C i ) + rkc( i ) 
r k d ( 1 ) = d t * ( v s / 1 1 ) 
rkd(3) s =dt*(-yyyyt4)/l£+vs/l£) 
t k d ( 4 ) ~ d t * ( y y ( 3 ) + y y y y ( 5 ) ) / c2 

rkd(5)*-dt*( (yyyy(4)+yyyy(6) )/lo+yyyy (5)*r/lo) 

rkd(6)=dt+yyyy(5>/cQ 

do 88 i = F6 

8S y ( i ) * y ( i ) + ( 1 .OdO/6.QdO)*(rka(i>+2.0dO*rkb(i)+e.OdO*rkc{l>+rkd<i) 

t = T"f-dt 

• i f ( y ( 4 ) . It . 0 . OdO ) then 
. y ( 4 ) * Q . Od 0 
endi f 

, £ s l = y ( l )-y C5) 

ss2=y (3)+*y (5) 
y(£)*0.0d0 

- 1 f ( m , ne . 1 ) go to 33 
i f ( n * e q * 1 0 ) then 
n = 0 

write(21,*)T,y(l> 
wr i x * < 22 , * > t , y ( 2 ) 
write(23,*)t,y(3) 

'HT it«(£4, 5 « i )tiy(4) 

. wri t e ( ZS t * > t , 4 y < 5 ) 

wri t e ( 26 / ^ ) t /y ( 6 ) 

w r i t e ( 2 7 . ^)t , ss! 
write(S^ < *)t,ss£ 
endif 



go to 300 
ON 


n-nT l 

33 i f ( i . 1 1 . ( t I + t 2 ft 3 ) ) 
c FOR CASE OF Si - ON A S2 - 


y ( 2) = 0 . OdO 

y ( 4 ) = o . o d o 

4 0 u r K a { t ) = d t * ( v s / 1 1 ) 

rka(3)=dt* (vs/12) 
rko(5)--dT j M y( 6)/lo*iy ( 5 ) * r / 1 o ) 
r k a ( S ) - d c : ^ ( y(5)/co) 
y y ( 1 ) = y ( 1 ) t-r k a ( t ) /£ . 0 dO 
yy ( 3)*y ( 3 ) +rka ( 3 )/2 . OdO 
yy ( 5)®y ( 5 ) +rka ( 5 )/2 . OdO 
yyC6)»y(6)+rka<6)/2. OdO 

r k b ( 1 ) * d t * ( v s / 1 1 ) 
rfcb ( 3 ) *dt ♦ ( vs/1 2 ) 

Pkb (5 ) »-dt * ( y y ( 6 ) / lo + yy (S ) +r/ lo ) 

rkb(6)«dt*( y y ( 5 )/co ) 

yyy(l )=yy( 1 ) -trkb ( 1 )/£ . OdO 

yy y ( 3 1 *=y y ( 3 ) +r kb ( 3 )/2 . OdO 

yyy(5)=yy(5)+rkb(5)/2. OdO 

yyy CS ) -y y ( 6 ) *rkb ( 6 )/2 . OdO 

rkc( 1 >=*dt* (vs/1 t ) 

.rkc(3)*dt* ( vs/12) 

rkc (5) =-dt * t yyy ( 6 )/i<H-y yy ( S) *r/l o ) 

rkc(6)=dt*( y y y ( S ) / co ) 

yy yy ( 1 ) ~y y y ( 1 ) +rkc ( 1 ) 

yy y y (3 )«y y y ( 3 ) +rkc ( 3 ) 

yyyy(5)»yyy(5)+rkc(5) 

yyy y (6 ) a =yyy(S)+rkc(6 ) 

r k a ( i) =dt * < vs/ 1 1 ) 

rkd ( 3 ) -d t ** ( vs / 1 2 ) 

rkd(5)=-dT*Cyyyy(6)/lo+yyyy(5)*r/lo) 

rkd(6)=dr^‘(yyyy(S)/co) 

y ( ( ) =y ( 1 ) t ( 1 . 0d0/6.0d0)*(rka(1 > +2 . OdO * r kb ( t ) +2 . OdO^rkc < 1 ) +r kd ( t ) ) 
y ( 3 ) = y ( 3 ) + ( i .0d0/6.0d0)*(rka(3)+-2.0d0*rkb(3)+2.0d0*rkc(3)*rkd(3) ) 
y ( 3 ) = y ( 5 ) f ( 1 . OdO/6. OdO ) * ( r k a ( 5 ) + £ . 0 d 0 * r k b ( 5 ) +2 . OdO*rkc(5)+rkd(5) ) 
y ( 6 ) - y ( 6 > + ( t . OdO/6. 0d0)*(rka(6>+£. 0d0*rkb(6)+£. OdO*rkc(6)+rfcd(6) ) 
T — T + d T 

ss 1 =y ( 1 ) -y ( S ) 

*. ss£-y(3)+y(5) 

1 f ( m . n e . 1 ) go to £2 
i f ( n . e q . 1 0 ) Then 
• . n~0 

write (21 , * ) t , y < I ) 

■wr i t e < 22 > *) t , y ( 2 ) 
writ e ( 23 , * ) t , y ( 3 ) 
write(£4,*)T./ y ( 4 ) 
wr it e ( £5 , * ) t , y ( S ) 
wr i t e ( 26 , * ) t , y < 6 > 
wr i t e ( 27 , * ) r , ssl 
writ e(£8j * ) t , ssH 
endif 
n = n+i 

2 i f ( t J1 e . ( 1 1 + 1 2 + 1 3+t 4 ) ) go to 400 

do 12 i a l,6 
yt2(i)*y(i) 

. err ( i ) =abs ( y 1 2 ( i )— y t 1 ( i ) ) 

]2 ' continue 

i f ( err ( 1 ) .gt. 0 . 1 d-8 ) go to 1000 
if (err (2) . g t . 0 . 1 d~8 ) go to 1000 

1 f ( err ( 3 ) . gt . 0 , 1 d-8 ) go to 1000 

if(err(4) .gt. 0.. Id-8) go to 1000 
i f ( err( 5 ) .gt, 0.}d~8) goto 1000 
if ( err ( 6 ) .gc. 0 . 1 d-8 ) go to 1000 



i f ( m . 1 1 . 2 } go to 1000 

clt>s?(8f ) ' 

c J o 5 e ( c 2 ) 

c3 ose( £3 ) 

cltt?e(c<!) 

cl 0 5S- (25 ) 

C 1 050 ( 2 6 ) 

*c J t> s o ( 2 7 ) 
cl os o(2P.) 

5 1. 0 p 
■? nd 



